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ABSTRACT 

We present the results of the first major systematic submillimetre survey of radio 
galaxies spanning the redshift range 1 < z < 5. The primary aim of this work is to 
elucidate the star-formation history of this sub-class of elliptical galaxies by tracing the 
cosmological evolution of dust mass. Using SCUBA on the JCMT we have obtained 
850-/im photometry of 47 radio galaxies to a consistent rms depth of 1 mJy, and have 
detected dust emission in 14 cases. The radio galaxy targets have been selected from 
a series of low- frequency radio surveys of increasing depth (3CRR, 6CE, etc), in order 
to allow us to separate the effects of increasing redshift and increasing radio power on 
submillimetre luminosity. Although the dynamic range of our study is inevitably small, 
we find clear evidence that the typical submillimetre luminosity (and hence dust mass) 
of a powerful radio galaxy is a strongly increasing function of redshift; the detection 
rate rises from ~15 per cent at z < 2.5 to >75 per cent at z > 2.5, and the average 
submillimetre luminosity rises at a rate oc (1 + z) 3 out to z ~ 4. Moreover our extensive 
sample allows us to argue that this behaviour is not driven by underlying correlations 
with other radio galaxy properties such as radio power, radio spectral index, or radio 
source size/age. Although radio selection may introduce other more subtle biases, 
the redshift distribution of our detected objects is in fact consistent with the most 
recent estimates of the redshift distribution of comparably bright submillimetre sources 
discovered in blank field surveys. The evolution of submillimetre luminosity found here 
for radio galaxies may thus be representative of massive ellipticals in general. 

Key words: galaxies: formation - galaxies: elliptical and lenticular, cD - dust, ex- 
tinction - radio continuum: galaxies - galaxies: active - stars: formation 



1 INTRODUCTION 

Although large numbers of 'normal' galaxies have now been 
discovered out to z ~ 5 (e.g. Steidel et al. 1999), radio galax- 
ies continue to offer the best opportunity to study exam- 
ples of massive elliptical galaxies (or the progenitors thereof) 
back to comparably early cosmic epochs. The reason for this 
is that a powerful radio source requires a massive black hole, 
and it is relatively certain that nowadays all such massive 
black holes reside in massive ellipticals (Magorrian et al 



mail: e. archibald@j ach. hawaii . edu 



1998; McLure et al. 199E). Thus radio selection offers a rela- 
tively efficient way of studying the properties of a sub-set of 
massive ellipticals as a function of cosmic epoch. Moreover 
this subset may well be representative of massive ellipticals 
in general, especially since it can be argued that a substan- 
tial fraction of all present-day ellipticals brighter than 2L* 
must have been active at z ~ 2.5 (Dunlop et al. 2000). 

At low redshift, radio g alaxies are domin ated by well- 
evolved stellar popul ations (Nolan ct al. 20001) , have rather 
low dust/gas masses (Knapp & Patten 1991), and lie in the 
same region of the fundamental plane as normal inactive el- 
lipticals (McLure et al. 1999; Dunlop et al. 2000). Moreover 
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current evidence suggests t hat they have evolved o nly pas- 
sively since at least z ~ 1 ( Lilly fc Longair 1984| ; McLun 
& Dun op 2000). This, coupled with the rela tively old age s 
derived for a few radio galaxies at z ~ 1.5 (Dunlop 1999), 
points towards a high redshift of formation, z > 3, for the 
bulk of their stellar populations. This means that both the 
star-formation rate and gas mass in these galaxies should be 
a strongly increasing function of redshift as one approaches 
their primary formation epoch(s). 

As has been argued by many authors, a massive star- 
burst at high redshift is expected to produce rapid chemical 
enrichment and to be largely enshrouded in dust. Conse- 
quently , submillimetre luminosity should be a good indi- 
cator of the evolutionary state of an elliptical galaxy, being 
expected to peak roughly half-way thro ugh the production 



of a galaxy from primordial material (Eales & Edmunds 
1996; Erayer & Brown 1997), or even earlier (i.e. as soon 
as the gas receives sufficient heating) given pre-existing en- 
richment (Hughes, Dunlop & Rawlings 1997). Indeed sub- 
millimetre luminosity, viewed as a tracer of gas mass, is ar- 
guably the best way to assess the evolutionary status of a 
massive galaxy. There are two reasons for this. Firstly, the 
sensitivity and bandwidth limitations of present-day instru- 
ments makes detecting molecular line emission from high- 
redshift objects extremely difficult. Secondly, in the case of 
radio galaxies, optical-UV measures may be confused by the 
direct or indirect effects of AGN activity. 

The detectability of high-redshift radio galaxies at sub- 
millimetre wavelengths was first demonstrated by Dunlop et 
al. (1994), when they detected 4C41.17, at the time the most 
distant known galaxy at z = 3.8. However the observation, 
made with the single element bolometer detector UKT14 
on the James Clerk Maxwell Telescope (JCMT), required 
4 hours of integration in exceptional weather conditions. In 
fact, the sensitivity of UKT14 permitted only the most ex- 
treme objects to be detected at high redshift, and it was 
often a struggle to detect even those. 

The advent of the Submillimetre Common-User 
Bolometer Array (SCUBA) on the JCMT offered the first 
real opportunity to rectify this situation. Although photo- 
metric observations do not fully exploit the multiplex advan- 
tage offered by an array camera, the individual bolometers 
in the SCUBA array offered almost an order of magnitude 
improvement in sensitivity over UKT14. This made it fea- 
sible to consider undertaking the first major submillimetre 
study of radio galaxies spanning a wide range of redshifts, 
and it is the results of the first major SCUBA survey of 
radio galaxies which we report here. 

The layout of the paper is as follows: In Section |^ we 
give a more detailed overview of pre-existing submillime- 
tre observations of radio galaxies, and explain the motiva- 
tion for observing a sample of galaxies compiled from flux- 
limited radio surveys of increasing depth. The resulting sam- 
ple is then described and summarized, before the results of 
our new submillimetre observations are presented in Sec- 
tion ^. Section ^ then gives details of the radio properties of 
each galaxy, and explains how the total and, where possible, 
core radio spectrum has been extrapolated to submillimetre 
wavelengths to estimate (or at least constrain) the potential 
level of non-thermal contamination at 850 /^m. The cover- 
age of the radio-luminosity:redshift plane provided by our 
observed sample is presented in Section ^, and then in Sec- 



tion we calculate the rest-frame 850-/im luminosities or 
upper limits for all the observed galaxies. In Section |^ we 
present a detailed statistical exploration of the evidence for 
genuine cosmological evolution of Lg5o Mm in our sample. Fi- 
nally in Section H we conclude by considering our results in 
the context of the recently published blank- field submillime- 
tre surveys. 

We have deliberately confined this paper to a determi- 
nation and discussion of the relative behaviour of submil- 
limetre luminosity in our sample, and have postponed our 
analysis and interpretation of more model-dependent prop- 
erties, such as inferred gas mass and galaxy age, to a subse- 
quent paper. 



1.1 Conventions 

For clarity, we summarize here a number of conventions 
which have been adopted throughout this paper. 

(i) For each galaxy in the sample, information has been 
collated from several references. Each reference has been 
given a code; for example, ER93 corresponds to Eales S.A., 
Rawlings S., 1993, ApJ, 411, 67. These codes are included 
in the bibliography. 

(ii) Throughout the paper, upper limits are calculated us- 
ing the following prescription: Consider an observation with 
signal S and standard error e. If the signal is positive, the 
n-a upper limit is S+(n x e). If the signal is negative, the 
n-a upper limit is n x e. It cannot be guaranteed that the 
upper limits taken from other papers were calculated in this 
manner. 

(iii) For a power-law spectrum, the spectral index, a, is 
defined as S v cx v~ a , where S„ is the flux density at fre- 
quency v. Thus, a radio synchrotron spectrum has a posi- 
tive spectral index, and the Rayleigh- Jeans tail of a thermal 
spectrum has a negative spectral index. 

(iv) The cosmological constant, A, is assumed to be zero 
throughout the paper. 



PROJECT BACKGROUND AND SAMPLE 
SELECTION 



2.1 



Pre-existing millimetre / submillimetre 
detections of high-redshift radio galaxies 



This work was carried out using SCUBA (Holland et al 



199S), the submillimetre bolometer array mounted on the 
JCMT. When the project was planned, SCUBA had not 
yet been commissioned, and only three high-redshift ra- 
dio galaxies had been unambiguously detected at mil- 
limetre/submillimetre wavelengths: 6C0902+34 (z=3.395), 
4C41.17 (z=3.8), and 8C1435+635 (z=4.25). While the mil- 
limetre observations of 6C0902+34 appeared to b e domi- 
nated by no n-thermal emission from the radio core ( Dowries 
et al. 1996| ), the observations of 4C41.17 and 8C1435+635 



were striking in their similarity, with high emission levels 
attributed to the presence of large reservoirs of dust. The 
observed flux densities are detailed in Table ^| 

4C14.17 and 8C1435+635 are both ultra-steep spec- 
trum (a > 1.0) radio galaxies, and both held the title, at 
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Figure 1. Radio luminosity-redshift plane as defined by radio surveys with successively deeper flux-density limits. Open circles denote 
galaxies from the 3C/3CRR survey and triangles denote galaxies from the 6CE survey. Stars indicate additional galaxies targeted as 
part of our SCUBA survey from the 4C*, 7CRS, 8C*, MIT-Green Bank MG*, and Texas TX* surveys. 4C41.17, 8C1435+635, and 
6C0902+34 are each indicated by a large cross surrounded by a circle. As can be seen from this figure, they are some of the most distant, 
radio-luminous objects known. Note, 53W069 is too faint to appear on this version (and all subsequent versions) of the P-z plane. H = 50 
kms _1 Mpc _1 , n a = 1.0. 



the time of their respective discovery, of being the most dis- 
tant radio galaxy known. In addition, as can be seen in Fig- 
ure |l|, they are two of the most-luminous (at 151 MHz) radio 
sources known. 

4C41.17 was d e tected at both 800 u m and 1300 fim 
( |Dunlop et al. 1994 |Chini fc Kriigel 1994 ), with the contri- 
bution from the non-thermal radio spectrum being less than 
a few percent at these wavelengths. Further support for a 
thermal emission mechanism came from the submillimetre 
spectral index: «?3oo ~ —4 (for self-absorbed synchrotron 
emission, a > —2.5). 



8C1 435+635 was first detected at 1250 /mi (Ivison 
199S|)~|Vs for 4C41.17, the contribution to the flux density 
from radio synchrotron emission is ~ 1 per cent. In addi- 
tion, the flux density is roughly the same as that detected 
for 4C41.17 at 1300 /im. 

Gi ven its similarity to 4C41.17, and its tantalizing de- 



suming H = 50 kms -1 Mpc and Q = 1). However, higher 
temperatures, and correspondingly lower dust masses, are 
not ruled out by the data. The far-infrared luminosity im- 
plied by the preferred fit is Lpir, ~ 1O 13 L.0. 

The amount of gas available for future star forma- 
tion can be calculated from the dust mass if a gas-to- 
dust ratio is assumed. This assumption is notoriously un- 
certain; however, robust upper and lower limits on the 
amount of gas present can be calculated if the highest and 
lowest reasonable values for the ratio are considered. For 
8C1435+635, such an analysis restricts the gas mass to 
4 x 10 10 < M gas < 1.2 x 10 12 M . In addition, if the dust 
is heated solely by star formation, the far-infrared luminos- 
ity indicates a star-formation rate of several thousand solar 
masses per year. 

The dust masses responsible for the detectio ns of 
4C41.17 and 8C1435+635 a re in excess of 10 s M 



tection at 1250 /im, 8C 1435+635 was a prime target for lop et al. 1994; Ivison 1995). This is at least an order or 



(Dt 



SCUBA. When this project was given its initial allocation 
of time, 8C1435+635 was the first object attempted, and 
we were able to study its submillimetre SED in unprece- 
dented detail. Firm detections were made at 850 (jm and 
450 /im, confirming that the submillimetre spectrum rises 
well above the declining radio emission. Observations at 
350 /im, 750 /im, and 175 /im resulted in upper limits to 
the continuum flux density; the last observation was made 
with the Infrared Space Observatory (ISO). In addition, a 
sensitive upper limit, 3<r < 5 x 10 10 Kkms _1 pc 2 was ob- 
tained for the CO(4— 3) line luminosity using the IRAM 30- 
m t elescop e. These data have been published in Ivison et 
al. ( 1998a| ); the observed flux densities (or limits) between 
A = 175/im and A = 1250/im are listed in Table |^. 

Isothermal fits to the submillimetre spectrum suggest 
the majority of the emission comes from 2 x 10 s Mq of dust 
with a temperature 40 ± 5 K and emissivity index j3 = 2 (as- 



magnitude higher than what is observed for radio g alaxies 
in the low-redshift Universe (Knapp & Patten 1991), indi- 



cating that, at z — 4, a significant amount of gas is yet to 
be converted into stars, suggesting th at the host galaxies are 
not fully formed (Hughes et al. 1997). 



2.2 Motivation for a complete sample of radio 
galaxies 

The distant radio galaxies initially detected at submillime- 
tre wavelengths are some of the most distant, radio-luminous 
objects known (Figure [j]) . On their own, they cannot be used 
to determine the evolutionary status of massive ellipticals at 
high redshift - their extreme submillimetre properties could 
be associated with either their redshift or the extreme na- 
ture of the radio source. There is a clear need for a proper 
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sample of radio galaxies spanning a range of radio luminosi- 
ties and redshifts. If the radio luminosity-redshift plane (the 
P-z plane) is properly sampled, the effects of radio luminos- 
ity and cosmological evolution can be disentangled. The goal 
of this project was to study such a sample using SCUBA. 

The galaxies were chosen based on their radio lumi- 
nosity and their redshift. Before doing so, it was necessary 
to determine the frequency at which the radio luminosity 
should be calculated. 

The radio luminosity of a classical double radio source is 
controlled by three key physical parameters: the bulk kinetic 
power in the jets, the age of the source (measured from the 
jet-triggering event) as it is observed on our light cone, and 
the environmental density (e.g. Blundell, Rawlings & Willott 
1999). Provided the rest-frame frequency is high enough so 
that it exceeds the synchrotron self-absorption frequency, 
and low enough that the radiative timescale of electrons in 
the lobes is not so short that the lobes are completely ex- 
tinguished, the selection effects inherent in any flux-limited 
sample lead to a reasonably close mapping between radio 
luminosity and jet power (e.g. Willott et al. 1999). These 
considerations render 151 MHz the rest-frame frequency of 
choice if one wishes to use radio luminosity to track the 
underlying jet power. 

In any flux-density limited sample, there is a tight cor- 
relation between luminosity and redshift: for large distances, 
only the brightest objects will be observable - the faint ob- 
jects will fall below the detection limit, and for small dis- 
tances there are few bright sources in the small available 
cosmological volumes. Suppose galaxies were selected from 
a flux-limited radio survey and observed with SCUBA, with 
larger submillimetre flux densities measured for the high- 
redshift galaxies than for the low-redshift galaxies. It would 
be impossible to determine whether this was genuine cos- 
mological evolution, or if the higher-redshift galaxies were 
brighter in the submillimetre because, for example, they are 
associated with objects with more powerful jets, perhaps 
because they are associated with more massive black holes 
and/or dark matter haloes. This problem can be overcome 
by using several radio surveys with different flux density 
limits to select our sample. The details of the redshift sur- 
veys of radio sources from which the sample was selected are 
described in Section 2.3. The region of the 151 MHz radio 



luminosity versus redshift plane covered by these surveys is 
shown in Figure |l]. 



2.3 The sample 

Initially, a sample of radio galaxies was compiled from the 
3C (Bennett 1962; Laing, Riley & Longair 1983 describe 
the doubly-revised version of the 3C sample, 3CRR) and 
6CE (Eales 1985; Rawlings, Eales & Lacy 2000) radio sur- 
veys, for which complete redshift information exists. The 
galaxies were chosen to span a range of redshifts, 0.2 < 
z < 4.5, but to lie in a narrow band of radio luminosity, 
10 27 WHz-'sr' 1 < Pisimhz < 10 29 WHz-'sr" 1 . This selec- 
tion was made to minimise any radio luminosity bias and to 
determine if the submillimetre emission from radio galaxies 
truly evolves with redshift. Note that the top two decades in 
radio luminosity were chosen, allowing galaxies to be stud- 
ied out to very high redshift. If there is a positive correlation 
between submillimetre emission and radio luminosity, per- 



haps via galaxy mass, this choice should also maximise the 
chance of making detections. 

In order to avoid strong radio synchrotron emission 
at submillimetre wavelengths, a further restriction was im- 
posed on the selected sources: they were required to have 
steep radio spectra, with a spectral index a > 0.5 for 
v ~ 1 GHz. This selection criterion was aided by the 
fact that the radio galaxies were selected to be very lu- 
minous at 151 MHz - the most intrinsically luminous ob- 
jects in low-frequency surveys tend to be steep-spectrum 
sources, whereas high-frequency surveys are biased towards 
flat-spectrum sources. 

The 6CE galaxies have a declination limit 5 < 40°. In 
order to minimise airmass, the criterion S < 40° was also 
imposed on the 3CRR catalogue sources. 

As the project progressed, it became apparent that 
sources selected from the 3CRR and 6CE samples alone 
could not produce an even spread in radio luminosity within 
the chosen luminosity band. Additional sources, often with 
declinations S > 40°, were added to the sample from: 

(i) two flux-density limited radio surveys with complete, 
or near-complete, redshift information: the 7C Redshift Sur- 
vey (7CRS) - Blundell et al. (2000); Willott et al. (2000); 
and the Leiden-Berkeley Deep Survey (LBDS) - Waddington 
et al. (2000) and references therein. 

(ii) several 'filtered surveys' for which a filter, spectral 
index and/or angular size for example, has been applied to 
identify high redshift galaxies in a radio survey. The names 
and references for these filtered surveys are as follows: 4C* - 
Chambers et al. (1996); 6C* - Blundell et al. (1998), Jarvis 
et al. in preparation; 8C* - Lacy (1992); TX* (Texas) - van 
Ojik et al. (1996); and MG* (MIT-Green Bank) - Stern et 
al. (1999), Spinrad private communication. 

Names, positions, and redshifts for the full sample ob- 
served with SCUBA are given in Table |l| In total, 47 radio 
galaxies were observed (27 from our original 3C/3CRR/6CE 
list). Note, for 3C356 two unresolved radio 'cores' have been 
detected within 5" of each other; they each have a galaxy at 
z=1.079 associated with them. It is not clear which galaxy is 
the host of the radio source; for further discussion of the mer- 
i ts of the two positions, refer to Best, Longair, & Rottgering 
( 1997 ) . SCUBA was pointed midway between the two identi- 
fications, and will have been offset 2" from the true position 
of the host galaxy; the flux lost should be negligible, < 5 
per cent at 850 /im. 



3 SUBMILLIMETRE OBSERVATIONS 

In this section, the results of all the submillimetre observa- 
tions of the 47 galaxies in the sample are presented. 



3.1 The SCUBA survey 

SCUBA contains two bolometer arrays cooled to ~ 75 mK to 
achieve sky background-limited performance. Each array is 
arranged in a close-packed hexagon with a 2.3' instantaneous 
field of view. The diffraction-limited beamsizes delivered by 
the JCMT to SCUBA are 14.7" at 850 /im and 7.5" at 
450 /im; the bolometer feedhorns are sized for optimal cou- 
pling to the beams. Consequently, the long-wavelength (LW) 
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Table 1. Positions and redshifts of the radio galaxy sample observed with SCUBA. The parent samples, described in Section \>. ] . arc 
listed in Column 5. Column 6 (Pos. ID) indicates how the position was measured: 'c' - a radio core ID exists, 'mc' - a radio ID exists 
for a marginal radio core, 'h' - the radio ID is the mid-point of the hotspot peaks, 'CSS' - same as 'h' but for an extremely compact 
radio source, 'O' - an optical ID exits, 'IR' - an infrared ID exists, 'HST' - a Hubble Space Telescope ID exists. The position and redshift 
references are given in column 8, with references containing redshift information in bold. Notes: § 6C0930+38 is also referred to as 
6C0929+38 in the literature. In the pre-release version of the 6C survey, astrometry indicated an RA of 09h 30m. However, in the final 
release of the survey, the position of the radio source appears as 09 29 59.7 +38 55 2. * The redshifts of 6C0919+38 and 6C1159+36 are 
not yet confirmed. For 6C0919+38, the estimate may not be too far off; for 6C1159+36, the redshift is based on a highly tentative Ly a 
emission line which is no longer considered secure (Rawlings et al. 2000). t For 3C356 the telescope was pointed at the midpoint of the 
two possible identifications (BLR97). * The redshift of MG1744+18 is unpublished. A reference to a paper which quotes the unpublished 
value is included. 



Common Name 


IAU Name 


RA (B1950.0) 


Dec (B1950.0) 


Parent 


Pos. 


z 


References. 




(B1950.0) 


(h m s) 


(° ' ") 


Sample 


ID 






6C0032+412 


0032+412 


00 32 10.73 


+41 15 00.2 


6C* 


c 


3.66 


BRE98, Jprep 


6C0140+326 


0140+326 


01 40 51.53 


+32 38 45.8 


6C* 


h 


4.41 


BRE98, RLB96 


4C60.07 


0508+604 


05 08 26.12 


+60 27 17.0 


4C* 


mc, O 


3.788 


CMvB96, Rv097 


4C41.17 


0647+415 


06 47 20.57 


+41 34 03.9 


4C* 


c, O 


3.792 


COH94, Rv097 


6C0820+36 


0820+367 


08 20 33.96 


+36 42 28.9 


6CE 


c, IR 


1.86 


ERLG97, LLA95, RELOO 


5C7.269 


0825+256 


08 25 39.48 


+25 38 26.5 


7CRS 


IR 


2.218 


ER96,W7CRS, BRR00 


6C0901+35 


0901+358 


09 01 25.02 


+35 51 01.8 


6CE 


IR 


1.904 


ERLG97, REW90, RELOO 


6C0902+34 


0902+343 


09 02 24.77 


+34 19 57.8 


6CE 


c, O, IR 


3.395 


COH94, L88, RELOO 


6C0905+39 


0905+399 


09 05 04.95 


+39 55 34.9 


6CE 


c, IR 


1.882 


LGE95, RELOO 


3C217 


0905+380 


09 05 41.34 


+38 00 29.9 


3CRR 


mc, HST 


0.8975 


BLR97, Lpc, RS97, SD84b 


6C0919+38 


0919+381 


09 19 07.99 


+38 06 52.5 


6CE 


IK 


1.65* 


ER96, RELOO 


6C0930+38§ 


0930+389 


09 30 00.77 


+38 55 09.1 


6CE 


h 


2.395 


NAR92, ER96, RELOO 


3C239 


1008+467 


10 08 38.98 


+46 43 08.8 


3CRR 


mc, HST 


1.781 


Lpc, BLR97, MSvB95 


MG1016+058 


1016+058 


10 16 56.82 


+05 49 39.3 


MG* 


c, O, IR 


2.765 


DSD95 


3C241 


1019+222 


10 19 09.38 


+22 14 39.7 


3CRR 


CSS, HST 


1.617 


BLR97, SD84b 


8C1039+68 


1039+681 


10 39 07.75 


+68 06 06.6 


8C* 


o 


2.53 


Lthesis 


6C1113+34 


1113+349 


11 13 47.64 


+34 58 46.6 


6CE 


h, IR 


2.406 


LLA95, RELOO 


3C257 


1120+057 


11 20 34.55 


+05 46 46.0 


3C 


O 


2.474 


Spriv, vBS98 


3C265 


1142+318 


11 42 52.35 


+31 50 26.6 


3CRR 


c, HST 


0.8108 


FBP97, BLR97, SJ79 


3C266 


1143+500 


11 43 04.22 


+50 02 47.4 


3CRR 


mc, HST 


1.272 


BLR97, LPR92, SD84b 


3C267 


1147+130 


11 47 22.07 


+13 03 60.0 


3CRR 


c, HST 


1.144 


BLR97, Lpc, SD84b 


6C1159+36 


1159+368 


11 59 20.94 


+36 51 36.2 


6CE 


CSS 


3.2* 


NAR92, RELOO 


6C1204+37 


1204+371 


12 04 21.75 


+37 08 20.0 


6CE 


IR 


1.779 


ERLG97, RELOO 


6C1232+39 


1232+397 


12 32 39.12 


+39 42 09.4 


6CE 


c, IR 


3.221 


ERD93, NAR92, REW90, RELOO 


TX1243+036 


1243+036 


12 43 05.40 


+03 39 44.5 


TX* 


c, O 


3.57 


vOR96, RMC95, Rv097 


MG1248+11 


1248+113 


12 48 29.11 


+11 20 40.3 


MG* 


O 


2.322 


SDS99 


3C277.2 


1251+159 


12 51 03.85 


+15 58 47.1 


3CRR 


c, HST 


0.766 


Lpc, McC97, SDM85 


4C24.28 


1345+245 


13 45 54.66 


+24 30 44.7 


4C* 





2.879 


CMvB96, Rv097 


3C294 


1404+344 


14 04 34.06 


+34 25 40.0 


3CRR 


c 


1.786 


MS90 


8C1435+635 


1435+635 


14 35 27.50 
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3C437 


2145+151 


21 45 01.58 


+15 06 36.1 


3CRR 


HST 


1.48 


BLR97, MSvB95 


MG2305+03 


2305+033 


23 05 52.31 
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array, optimised for an observing wavelength of 850 /xm, 
contains 37 pixels. The short-wavelength (SW) array, opti- 
mised for operation at 450 (im, contains 91 pixels. A dichroic 
beamsplitter allows observations to be made with both ar- 
rays simultaneously. 

For point sources, SCUBA's photometry mode is rec- 
ommended, where the target is observed with a single pixel. 
During commissioning, it was discovered that the best pho- 
tometric accuracy could be obtained by averaging the source 
signal over a slightly larger area than the beam. This is 
achieved by 'jiggling' the secondary mirror such that the 
chosen bolometer samples a 3 x 3 grid with 2" spacing cen- 
tered on the source. The on-source integration time is 1 sec- 
ond per jiggle point; excluding overheads, the jiggle pattern 
takes 9 s to complete. While making the mini jiggle map, the 
telescope was chopped at 45" in azimuth at a frequency of 
7Hz. After the first 9-point jiggle, the telescope was nodded 
to the reference position, an event which occurred every 18 
seconds thereafter. 

Between June 1997 and March 1999, a total of 192 
hours was used to carry out the observations, all in excellent 
weather conditions. The atmospheric zenith opacities were 
at the very most 0.3 at 850 micron and 2.5 at 450 micron, al- 
though they were frequently much better than this. Skydips 
and observations of calibrators were obtained regularly. For 
each source, the data have been fiatfielded, despiked, and 
averaged over 18 seconds. The residual sky background that 
chopping and nodding are unable to remove has been sub- 
tracted by using the data taken by the off-source bolome- 
ters. The Kolmogorov-Smirnov (KS) test has been applied 
to the data to check for consistency, and the data have been 
calibrated using observations of Mars, Uranus and several 
secondary calibrators. 

The 850 /im and 450 /im flux densities measured for the 
radio galaxies in the sample are listed in Table ^. This table 
includes both 3a upper limits for sources whose signal-to- 
noise ratio (S/N) does not exceed 3.0 and 2a upper limits for 
sources whose S/N does not exceed 2.0. The reason for this is 
that while traditionally the detection of a given source might 
only be regarded as robust if it is greater than 3a, genuine 
2a detections should be taken seriously, especially given that 
the telescope was pointed at known galaxies instead of the 
submillimetre source being identified in a blank-field sur- 
vey. Furthermore, the adoption of a 2a detection level is 
necessary to implement the statistical analysis presented in 
Section [7.l| . It should be noted, however, that the results 
presented in this paper hold whether a 3a or a 2a detection 
threshold is adopted. 



observed with SCUBA (for 3C3 24, alt hough we find a lower 
flux density than Best et al. (1998b), the two results are 



equivalent within errors). There are two exceptions: 

(i) 53W002 - Ev en though 53W002 w as detected at 800- 
fj,m with UKT14 ( [Hughes et al. 1997[ ), we have failed to 
detect it with SCUBA. Assuming a typical submillimetre 
spectral index of —4.0, the 800-^im detection of 53W002 im- 
plies a flux density of at least 3.6 mjy at 850-^tm. Note the 
spectral index is expected to lie between the values of —3.0 
and —4.0 (e.g. Hildebrand 1983), and adopting a value of 
—4.0 will give the lowest estimate of 850-/xm flux. Thus, if 
the UKT14 data are trustworthy, we should have detected 
53W002 with SCUBA, even if the large error bar on the 
800- ^m measurement is taken into account. 

(ii) MG1016+058 - For MG1016+058, previous detec- 



tions e xisted at both 1300 /im and 800 fim (Cimatti et al 



1998b ). The minimum submillimetre spectral index consis- 
tent with these data points and their error bars is ~ —2.8. 
This can be used to estimate a lower limit on the flux den- 
sity expected at 850 /im. The 850-/im flux density should 
be well in excess of 5 mjy, easily detectable with SCUBA. 
An obvious objection to this analysis is that the 1300-/im 
flux density could be solely due to radio synchrotron emis- 
sion. This is unlikely; the radio spectrum is highly curved 
and dives down well below the 1300-/im flux density (Ap- 
pendix |X|). However, assume for a moment that the 1300- 
(im flux density was contaminated by radio emission. Again, 
adopting a typical submillimetre spectral index of -4.0 and 
applying it to the 800 /im detection still predicts an 850-/im 
flux density in excess of 5 mjy. 

These two discrepancies are perhaps not all that sur- 
prising. UKT14 was a single-element bolometer, relying 
solely on chopping and nodding for sky removal. However, 
a chop/nod set-up is unable to account for quick variations 
in the background sky, and there will be residual sky noise 
present in the data (Jenness, Lightfoot & Holland 1998). 
SCUBA is an array, and the off-sourc e bolometers can b e 
used to remove this r esidual sky noise (Jenness et al. 1998). 
Omont et al. (1996) investigated the difference between 
single- and multi-bolometer detectors, and found that the 
less-reliable sky-cancellation offered by single-bolometer sys- 
tems often resulted in fake detections of weak sources. Fur- 
thermore, they failed to detect three high-redshift quasars 
with a mult i-bolo meter array that Andreani, La Franca 
& Cristiani (1993) had apparently detected with a single 
bolometer. 



3.2 Radio galaxies with other 

millimetre / submillimetre observations 

Several attempts were made to observe the millime- 
tre/submillimetre spectrum of radio galaxies before SCUBA 
existed. The attempts, made with the IRAM 30-m telescope 
and UKT14 (SCUBA's predece ssor) , w ere largely unsuccess- 
ful. More recently, Best et al. ( |l998b|) also looked at 3C324 
with SCUBA and Benford et al. (119991) observed 4C41.17 at 



350 ^m with the CSO. All of these observations, including 
the full submillimetre spectrum of 8C1435+635, are listed 
in Table |[ 

In most cases, these data are consistent with what we 



4 CORRECTING RADIO GALAXY DUST 
EMISSION FOR POTENTIAL 
CONTAMINATION 

For the majority of galaxies in the sample which we have 
actually detected, a detection has only been obtained at 
one submillimetre wavelength (i.e. 850 /im). Without a reli- 
able estimate of the submillimetre spectral index, care needs 
to be taken to ensure the detected submillimetre emission 
is clearly in excess of the radio synchrotron spectrum, and 
can thus be safely attributed to dust. However, as will be- 
come clear, due to the fact that the sample was confined 
to steep-spectrum lobe-dominated radio galaxies at z > 1, 
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Table 2. Observed submillimetre flux densities (S„) and standard errors for the sample. 3cr upper limits are shown for sources whose 
S/N does not exceed 3.0; 2<r upper limits are shown for sources whose S/N does not exceed 2.0 
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for almost all objects, detectability with SCUBA at 850 
corresponds to a successful detection of dust. Radio spectra 
have been compiled for each source with the aim of extrap- 
olating the radio emission to submillimetre wavelengths and 
subtracting it off. This correction should give a good esti- 
mate of the flux density produced by thermal dust emission. 

It might appear that such a correction would be ex- 
cessive, or even completely unnecessary where it is known 
that the radio lobes lie outside the SCUBA beam. However, 
at high frequencies, > 8 GHz, experience with individual 



well-studied objects (e.g. 4C41.17) indicates that much of 
the emission can arise from structures coincident or close to 
the radio core (e.g. Carilli, Owen & Harris 1994). Therefore 
we have adopted the conservative but consistent approach 
of estimating submillimetre synchrotron contamination from 
the extrapolation of the total high-frequency radio emission 
in every case. In practice, for radio sources larger than the 
JCMT beam, this correction makes a significant impact on 
the estimated dust mass in only a very small n umb er of 
sources. These are discussed at the end of Section 4.1. 
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Table 3. All other millimetre and submillimetre observations of galaxies in the sample. All flux densities are quoted in mjy. The references 
are: 3C324 - BRB98; 53W002 - HDR97; 3C257 - HDR97; MG1016+058 - CFR98; 6C1232+39 - C&K94; 6C0902+34 - C&K94, DSS96, 
HDR97; TX1243+036 - CFR98; MG2141+192 - HDR97; 6C0032+412 - HDR97; 4C41.17 - BC099, C&K94, DHR94; 8C1435+635 - 
HDR97, 195, IDHA98 



Source 


S3000pm Sl300/jm 




S800/jm 


S750^m S45fj^m 




3C324 




3.65 ±1.17 




3(7 < 21 




53W002 






6.9 ± 2.3 






3C257 






3(7 < 11 






MG1016+058 


2.13 ± 0.47 




14.7 ±4.6 






6C1232+39 


3ct < 3.0 










6C0902+34 


4.2 ±0.6 3.1 ±0.6 




3(7 < 14 






TX1243+036 


3<r < 2.6 




3cr < 9.3 






MG2141+192 






3(7 < 11 






6C0032+412 






3(7 < 14 






4C41.17 


2.5 ±0.4 




17.4 ±3.1 


3(7 < 56 


37 ±9 


8C1435+643 


2.57 ±0.42 


7.77 ±0.76 


3cr < 13 


8.74 ±3.31 23.6 ±6.4 


3(7 < 87.0 3(7 < 40.1 



The radio spectra are shown in Appendix |X|, and include 
both integrated and core radio flux densities. Core flux den- 
sities are generally much fainter than the integrated radio 
emission at low frequencies, but they can also be much flat- 
ter. If the core has a flat spectral index, it may contribute 
significantly to the submillimetre flux density even if the 
integrated radio emission seems to fall well below this. 

Fits to the radio spectra are also shown in Appendix [X| 
These fits are described in Section [ll] (with full details in 
Appendix ^) and are subsequently used to estimate the non- 
thermal contamination at submillimetre wavelengths. De- 
tails on how the ind ividu al radio spectra were compiled are 
given in Archibald fll99S|) . 

Note the radio spectrum of 53W069 has not been pre- 
sented here. 53W069 is only 3.73 mjy at 1.4 GH z, and syn- 
chrotro n c ontamination will not be a problem (Windhorst 
et al. 1§8^). 



For 1.5 GHz < v < 5 GHz, the spectrum steepens to an 
index of a ~ 0.9. At high frequencies 5 GHz < v < 15 GHz, 
the spectrum steepens even further to an index of a ~ 1.0. 



4.1 Synchrotron contamination 

A simple synchrotron spectrum is a power law of the form 
S v oc v~ a . However, real radio spectra turn over at low fre- 
quencies owing to synchrotron self-absorption, or through 
running out of electrons at low energies. At high frequen- 
cies the spectra steepen as they age - the highest-energy 
electrons radiate away their energy the fastest. In log-space, 
such spectra appear curved. 

Consider, as an example, the radio spectrum of the 
hotspots of Cyg nus-A , described in detail by Muxlow, Pel- 
letier & Roland fll988j ). The spectrum is a power law at low 
frequencies 0.4 GHz < v < 1.5 GHz, with a spectral index 
a ~ 0.5, and a low-frequency turn-over at v ~ 0.2 GHz. At 
2GHz the spectrum steepens; for v > 2 GHz, the spectrum 
is cons istent with a power law of index a ~ 1.0. R obson 
et al. (1998) and Eales, Alexander & Duncan ( 198S| ) have 
found that this steep power law is an accurate description of 
the spectrum right out to submillimetre wavelengths. Thus 
in log-space, at high frequencies a straight-line fit would be 
accurate. However, when fitting to the whole spectrum, a 
model with some curvature would be required. 

Another good example is the study of the dom inant 
hotspot of 3C273, by Meisenheimer & Heavens ( 1986 ). For 
low frequencies, v < 1.5 GHz, the spectral index is a ~ 0.7. 



More recently, Murgia et al. (1999) published a study 
of compact steep-spectrum radio sources. The majority of 
objects in their sample exhibit a clear spectral steepening, 
occurring anywhere from a few hundred MHz to tens of GHz. 

In order to successfully correct for synchrotron contam- 
ination at submillimetre wavelengths, the high-frequency ra- 
dio spectrum needs to be well determined. Unfortunately, for 
most galaxies in the sample, the highest-frequency observa- 
tion is v ~ 20 GHz. A straight-line extrapolation (in log- 
space) from the highest frequency data point could be used 
to decontaminate the submillimetre flux densities. However, 
in the interval between 20 GHz and the SCUBA wavebands, 
the spectrum may steepen, and the straight-line extrapola- 
tion would over-correct the SCUBA flux density. 

Evidence to support this steepening comes from 3C324, 
MG1744±18, and MG1016±058 - the only sources in the 
sample with measured radio flux densities at v > 20 GHz. 
In each case the high-frequency observation indicates a 
significant steepening of the radio spectrum. 4C23.56 and 
8C1435±635, two of the higher-redshift galaxies in the sam- 
ple, even appear significantly curved for v < 20 GHz. 

Note that if a straight-line extrapolation of the high- 
frequency radio emission is correct, three of 850- /im detec- 
tions made with SCUBA appear to be seriously contami- 
nated with non-thermal synchrotron - 4C13.66, 3C470, and 
3C257. 

More high-frequency radio observations are needed to 
be certain what the radio spectrum of each galaxy in the 
sample does as it approaches the submillimetre waveband. 
For now, both possibilities will be considered: a linear fit, 
a straight-line extrapolation of the high-frequency data, as 
could be accurately applied to Cygnus-A; and a parabolic fit 
which reflects a degree of steepening, or curvature, at higher 
frequencies. The details of the two fits, together with plots 
of the SEDs, are given in Appendix 

It would be surprising if the radio spectra flattened off 
at high frequencies, as this only occurs if the source has an 
extremely bright, flat radio core (except for 6C0902±34, the 
radio gal axies discussed here all have relatively faint cores 
- Section |4.2[). Thus, the linear fit is a strong upper limit 
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on the synchrotron emission at submillimetre wavelengths, 
especially since a significant fraction of the emission may 
lie outside the JCMT beam. Deviations from this standard 
power-law take the form of a steepening at high frequen- 
cies - the parabolic fit mimics this curvature and is a good 
lower limit on the amount of synchrotron emission at submil- 
limetre wavelengths. Given the lack of high-frequency radio 
observations, it is often difficult to determine which is the 
better fit. The best estimate of the synchrotron contamina- 
tion has been taken to be the midpoint of these two models, 
with the error in the estimate being the difference between 
the midpoint and either the upper or lower limit. 

These estimates of synchrotron emission were thus sub- 
tracted from the measured SCUBA flux densities, with the 
errors being added in quadrature (as prescribed by the prop- 
agation of errors) . However, there were four special cases for 
which particular care had to be taken: 

(i) 3C324, MG1744+18, MG1016+058, and 4C23.56 all 
have high-frequency data that the linear fit has trouble fit- 
ting to. For these galaxies, the parabolic fit looks to be the 
more realistic extrapolation of the radio spectra, and was 
used to decontaminate the SCUBA flux densities. 

(ii) For the 14 radio sources whose lobes lie outside the 
JCMT beam, we may have over-corrected for synchrotron 
contamination. However, for all of these but one, the source 
was either undetected by SCUBA, or the applied correction 
was insignificantly small. For 3C470, the applied correction 
turned a detection into a non-detection, reducing the 850- 
/j,m flux from 5.6 mjy to 3.1 mjy. However, as mentioned 
previously, our method of correction may be reasonable for 
sources that lie outside the JCMT beam. Furthermore, if 
the correction is not applied to 3C470, the results presented 
in the paper are unaffected. 

(iii) For some galaxies SCUBA measured a negative sig- 
nal. They have not been corrected for synchrotron contami- 
nation; performing the correction would have unnecessarily 
increased the noise of the measurement and inappropriately 
made the signal even more negative. 

(iv) For some of the undetected galaxies, the estimated 
synchrotron contamination is larger than the signal mea- 
sured by SCUBA, resulting in a negative submillimetre flux 
density if it is corrected for. Claiming a 'negative emission' 
from dust is not physically meaningful. It is better to set the 
SCUBA flux densities equal to zero: it appears that there is 
no thermal emission from dust at these wavelengths. This 
leaves the question of how to handle the errors. It seems 
best to leave them as they are - as they have been measured 
by SCUBA. The statement then being made is that, as far 
as it can be ascertained, all of the measured signal may be 
radio synchrotron emission, and any dust emission from the 
source is definitely less than 2x the measured noise (for a 
2<j upper limit). This seems a good, conservative upper limit 
- in each case the limit is larger than the original measured 
flux density. 



4.2 Radio cores 

When searching for thermal dust emission, the telescope is 
pointed at the optical/infrared ID of the galaxy. This posi- 
tion is coincident with the radio core (if one has been de- 
tected). A flat radio core could contribute at submillimetre 



wavelengths. A useful exercise is to estimate the upper limit 
on this contribution for each source - if it tends to be negligi- 
ble, the SCUBA flux densities will not need to be corrected. 

Only 8 of the galaxies in the sample have their radio 
core detected at more than one frequency: 3C267, 3C324, 
3C241, 4C23.56, 6C1232+39, 6C0902+34, TX1243+036, 
and 4C41.17. 6C0902+34 is a special case that is excluded 
from the following discussion; it will be treated separately 
at the end of this section. For the remaining 7 galaxies, the 
radio core spectral index has been estimated using the two 
highest-frequency core detections. 3C324 and 3C241 have 
relatively flat core spectra, with a ~ 0.35. The majority, 
however, have very steep spectra, ranging from a ~ 0.7 to 
a ~ 1.8. Note that no detected core appears either com- 
pletely flat or inverted (i.e. flux density increasing with 
frequency) between 1 GHz and 15 GHz. If these cores 
are extrapolated to submillimetre wavelengths, in all cases, 
save one, the core contribution is negligible: < 0.05mJy at 
850 fim. For 3C241, on the other hand, the core strength at 
850 (im is likely to be ~ lmjy. 

For several galaxies in the sample, the radio core has 
been detected at a single frequency. Adopting the flattest 
measured spectral index of a ~ 0.35, the corresponding core 
flux densities at 850 fim have been estimated. This is a pes- 
simistic estimate, many of these galaxies may have steeper 
core spectra, but it is a good indication of the worst-case 
scenario. In all cases but three, the predicted core flux den- 
sities are negligible, < 0.2mJy. Higher 850 fim core flux 
densities are expected for 3C322 and 3C265. However, for 
these galaxies a negative signal was measured with SCUBA, 
and the core brightness is irrelevant. For MG1016+058, if 
the spectrum is flat, a core flux density of ~ 0.8mJy is ex- 
pected at 850 fim. If, on the other hand, the spectral index 
is steeper, a ~ 0.75 for example, the core contamination is 
negligible. 

A lack of sources with core detections at more than one 
frequency potentially hinders this analysis. However, given 
the available information it seems unlikely that the cores 
are significant except for 3C241 and MG1016+058. The core 
contribution in the case of 3C241 has been corrected for, as 
the core spectrum is definitely flat. For MG1016+058, on the 
other hand, it is impossible to know whether the core is flat 
or steep. Given this uncertainty, the SCUBA observations 
have not been corrected for core contamination. Instead, it 
is simply noted that MG1016+058 may have a large core 
contribution at submillimetre wavelengths. 

It is mentioned in the literature that three galax- 
ies have completely flat (or possibly inverted) cores at 
v < 5 GHz: 3C294 (MS90), 6C0905+39 (LGE95), and 
6C0032+412 (BRE98; Jarvis et al., in preparation). These 
cores could contribute significantly to the observed S CUBA 
flux densities. However, Rudnick, Jones & Fiedler (1986) 
have studied flat cores in powerful radio galaxies. They 
found the cores often turned over or steepened at v > 5 GHz 
- a completely flat core at low frequencies does not mean 
that it will remain flat out to the submillimetre wavebands, 
indeed all evidence suggests it is very unlikely to do so. 

In addition, note that for 3C356, two radio core iden- 
tifications exist: one faint and steep, the other brighter and 
flat. Neither should significantly contribute at submillimetre 
wavelengths. 

6C0919+38 could have a very bright radio core, ~ 
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6 mjy at 5 GHz. However, it is not c lear w hether the bright 
feature detected by Naundorf et al. (1992) is the radio core 



the radio galaxy itself, thi s can not be regarded as a prob- 
lem. See also Blain et al. (1998) for an in-depth discussion 



or a knot in the radio jet. Given the uncertainty, the possi- 
ble radio core contamination is merely noted here and is not 
corrected for. 



4.2.1 6C0902+34 

6C0902+34 is a special case because of the brightness of its 
core, ~ lOmJy, between 1.5 and 15 GHz. In addition, two 
slightly conflicting observations at 8 GHz make it difficult 
to be certain how steep the core spectrum is. If the flatter 
spectral index is used (a ~ 0.3 ), the 3 mm and 1. 3 mm 



detections of Downes et al. ( 1996 ) and Chini & Kriigel ( 1994 ) 
look t o be c ompletely dominated by the radio core. Downes 
et al. (1996) came to the same conclusion, and predict the 
contribution of dust emission at 1.3 mm to be < 0.6 mjy. 
Even if the steeper spectral index is used (a ~ 1.0), the 
detections could easily be dominated by the integrated radio 
spectrum (although higher-frequency radio detections are 
needed to confirm this). 

When observed with SCUBA at 850 /im, a 3a upper 
limit of 5.83 mjy was measured. Combining this with Chini 
& Kriigel's detection of 3.1 ± 0.6 mjy at 1300 /im, the 
submillimetre spectral index is |a| < 2.0. This is inconsis- 
tent with thermal dust emission being present, for which 
|a| > 2.0. (Note that the submillimetre spectral indices are 
negative) . 

6C0902+34 appears completely dominated by radio 
synchrotron emission. It has thus been left out of the sam- 
ple. Hereafter, we discuss a sample of 46 radio galaxies which 
excludes 6C0902+34. 



of source confusion at SCUBA wavelengths. 



4-3.3 Confusion by sources in the off-beam 

The measured flux densities could be affected by chopping 
onto a nearby source. This is unlikely to be an impor- 
tant problem for this study, even though SCUBA maps of 
the fields of radio galaxies have revealed an over-density of 
sources in some cases (Ivison et al. 2000a; Ivison et al., in 
preparation). For the galaxies in our sample which have been 
mapped and shown to have companion sources (4C41.17 and 
8C1435+635), our chop did not land on the companions. 
Furthermore, our observations were made with an azimuthal 
chop, and the off-beam will have moved across a substantial 
angle (several degrees) of the sky as we tracked a source, di- 
luting the effect of off-beam contamination. There is also no 
a priori reason why this possible source of confusion would 
result in a redshift bias and affect the analysis presented in 
Section [?]. 



4.4 A final note on corrections: 6C0140+326 

As a final note on correcting the SCUBA flux densit ies, it ap- 
pears that 6 C0140+326 is gravitationally lensed (Rawlings 
et al. 1996a] ). Jowever, the predicted amplification factor is 
small (< 2). Given the uncertainties in estimating this fac- 
tor and the synchrotron contamination, a correction has not 
been applied. 



4.3 Submillimetre confusion 

4-3.1 Galactic cirrus confusion 

Extrapolating the galactic cirrus confusion that IRAS mea- 
sured at 100 /im to 850 /im at SCUBA resolution yields a 
cirrus confusion noise of << 1 mjybeam - in the direction 
of our targets. Thus, contamination by galactic cirrus should 
not be a problem for this study (see Hughes et al. 1997 and 
Helou & Beichman 1990 for detailed discussions). 



4-3.2 Confusion by extragalactic sources 



Hughes et al. ( 199S ) performed the deepest submillime- 



tre survey when they mapped the Hubble Deep Field with 
SCUBA. For 850 /im, confusion was found to be a prob- 
lem for sources weaker than 2 mjy. At the 2 mjy level, the 
source density of their map is ~ 1 source per 30 beams. 
In a large blank-field survey reaching 2 mjy, some detec- 
tions will inevitably result from confusion, particularly if 
sources are strongly clustered. However, the chances that 
such false confusion-induced detections should coincide in 
position with known massive objects at high redshift is 
clearly extremely small. Of course, given sufficient resolution 
our detected sources may break up into sub-components, as 
for example has the optical im age of 4C41.17 under HST 
resolution (Miley et al. 1992). However, since such sub- 
components are still extremely likely to be associated with 



4.5 Corrected flux densities 

The final 850-/im flux densities, corrected for synchrotron 
contamination, and in the case of 3C241 for radio core con- 
tamination, are presented in Table [|. Note, we do not in- 
clude the corrected 450-/im flux densities; at 450 fim the 
corrections are tiny compared with the errors. As mentioned 
earlier, we adopt a detection threshold of S/N> 2.0, with 
upper limits taken at the 2a level. Out of a sample of 46 ra- 
dio galaxies (excluding 6C0902+34), thermal emission from 
dust has been detected in 14 galaxies at 850 /im, and in 5 
of these galaxies at 450 /im. 

For 3C368, MG1016+058, 4C60.07, 4C41.17, and 
8C1435+635, an estimate of the 850-450 fim spectral index 
can be made. For thermal dust emission, the slope along the 
Rayleigh- Jeans tail is expected to be around 3.0-4.0. The 
spectral indices of the lower-z detected sources, 3C368 and 
MG1016+058, are consistent with this: |a||o| ~ 4. The in- 
dices of 4C60.07, 4C41.17 and 8C1435+635 are on the low 
side, |o4io| < 2. However, this is to be expected. They have 
large redshifts, and for thermal greybody emission at 30 K, 
the spectral turnover should be redshifted into or near the 
450 /im waveband. Using the 450 fim observation would un- 
derestimate the slope along the Rayleigh- Jeans tail in this 
situation. 
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Table 4. 850-/xm flux densities (S„) and errors corrected for ra- 
dio synchrotron contamination. In the case of 3C241, radio core 
contamination has also been taken into account. A '</' in column 
2 indicates that applying the correction has changed the flux den- 
sity by less than la. 2a upper limits are shown for sources whose 
S/N does not exceed 2.0. 



Source 



< la 



(mjy) 



850 Atm 

S/N 



2a limit 
(mjy) 
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5 THE RADIO LUMINOSITY- REDSHIFT 

PLANE FOR RADIO GALAXIES OBSERVED 
WITH SCUBA 

It is useful to plot radio luminosity against redshift (the 
P-z plane) highlighting the galaxies we have observed with 
SCUBA. This ensures that even luminosity coverage of the 
P-z plane has been achieved, and identifies the redshifts 
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Figure 2. Radio luminosity-redshift plane for Q a = 1.0 (upper 
plot) and f2 = 0.1 (lower plot), assuming H = 50 kms — 1 Mpc — 1 . 
The solid circles indicate radio galaxies for which thermal emis- 
sion from dust was detected at 850 /im with S/N> 2, and the 
open circles indicate galaxies for which no dust emission was de- 
tected. For the detections, the size of the solid circle represents 
the brightness of the galaxy at 850 fim. The 151-MHz radio lumi- 
nosities, Pi5iMHz, were estimated using the radio SEDs presented 
and referenced in Appendix Note, 53W069 is too faint to ap- 
pear on this version of the radio luminosity-redshift plane. The 
dashed-line boxes depict Subset A (upper plot) and Subset B 
(lower plot) which are defined in Section 7.1. 



at which submillimetre emission from radio galaxies is de- 
tected. 

As previously mentioned, the radio luminosities have 
been calculated at rest-frame 151 MHz to ensure that the 
radio luminosity tracks the power of the radio jets as ac- 
curately as possible. Simple physical models for the time 
development of the 151 MHz luminosity and the linear size 
(Diinear) of a double radio source suggest that over a large 
range of D; inear , the 151 MHz radio luminosity gives a rea- 
sonably accurate guide to the jet power Q of the vast major- 
i ty of objects detect ed in complete radio flux-limited samples 
(|Willott et al. 199§), 



Both the 151 MHz luminosities and the linear sizes for 
the sample have been calculated for two cosmologies: fi = 
1.0, and Q, = 0.1, where H =■ 50 kms _1 Mpc _1 . The values 
are presented in Table |^. Figure ^ shows the corresponding 
P-z plane for each cosmology. 
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A proper statistical analysis of these data will be pre- 
sented in Section [?], but at first glance, Figure || seems to 
indicate that submillimetre emission is more predominant in 
high-redshift radio galaxies than in low-redshift radio galax- 



6 SUBMILLIMETRE LUMINOSITIES AND 
DUST MASSES 

To calculate the rest-frame 850- /im luminosity (Lsso^m) and 
dust mass (Mdust) of each source, we adopt an optically 
thin isothermal dust emission template with (3 = 1.5 and 
Tdust = 40 K. This template is consistent with the studies 
attempting to constrain j3 and Td us t that have been pub- 
lished in the literature. For example, Benford et al. (199!:) 
have observed a sample of high-redshift radio-quiet quasars 
and radio galaxies at 350-/im. Assuming a critical wave- 
length of A = 125/xm, isothermal fits indicate /3 ~ 1.5 and 
Tdust ~ 50K. The study of IRAS galaxies by Dunne et al. 
1.2 and T dusrt 



( g000[ ) indicates /3 ~ 1.2 and T dus t ~ 36K (although the 
value /3 is thought to have been artificially lowered by the 
presence of a cold dust component) . Finally, for the handful 
of galaxies identified in the recent SCUBA surveys that have 
known redshifts (refer to Smail et al. 2000 for a review), the 
favoured value of Td us t is ~ 40 K (e.g. Ivison et al. 1998b, 
2000b; Barger, Cowie & Sanders 1999). 

It is currently impossible to determine the precise val- 
ues of f3 and Td us t, and even though we have assumed the 
most likely template, it could very well be wrong. This has 
important consequences for the analysis presented here. For 
greybody emission, the values of (3 and Td us t determine the 
shape and position of the spectral turnover, but have less 
of an effect on the Ray leigh- Jeans tail. As the redshift of 
the source increases, the spectral peak is brought into view. 
Thus the K-correction used to calculate the 850-/im lumi- 
nosity is more sensitive to changing the values of and Td us t 
at high redshift than it is at low redshift. If the template is 
wrong, the error in the estimations of Lsso^m will be larger 
if the source has a high redshift. To explore this, the analy- 
sis of evolutionary trends has been conducted using several 
/3-Tdust combinations, in addition to the template itself, to 
ensure the trends are real. 

Furthermore, applying a template to the entire sample 
assumes that the individual radio galaxies all have similar 
dust parameters. This is a reasonable assumption if the same 
mechanisms for creating and heating dust are expected in 
each source. If, however, the dust characteristics do vary 
from source to source, they would have to do so in a redshift- 
(or radio power-) dependent manner to seriously affect the 
evolutionary trends presented here. 

The choice of cosmology is also an important issue. As 
the relative properties of the sample are of most interest in 
the analysis presented here, the value of the Hubble constant 
is irrelevant. Changing the Hubble constant simply scales all 
the luminosities up or down by the same amount. However, 
changing the value of the density parameter affects the lumi- 
nosities in a redshift-dependent manner. In order to encom- 
pass all possibilities, we thus consider two extreme values: 
fio = 1.0 and D. = 0.1, with H D = 50 kms^Mpc -1 . In 



E 



0.01 




Section 7.2 a low-density Universe (fl = 0.1) with H = 67 
is also considered. The reason for this last choice is that the 



Figure 3. The 850 /jm flux density that would be observed for 
the starburst galaxy M82 (Hughes, Gear & Robson 1994) if it 
were placed at progressively higher redshifts, assuming the dust 
emission template adopted here (J3 = 1.5, T= 40K). The solid 
line represents an Q a = 1.0 Universe; the dashed line represents 
an Uo = 0.1 Universe (H = 50 km-'Mpc -1 ). 



age of a Universe with Q Q = 1.0, H = 50, is ~ 13 Gyr. A 
Universe with fl a — 0.1, H = 67 has the same age - hence 
the 'absolute values' of the luminosities in these high and 
low density cases can be realistically compared. Recall the 
cosmological constant, A, is assumed to be zero throughout. 
It is worth noting that the evidence from the Supernova Cos- 
mology Project is currently in favour of a cosmology with 
f2 M = 0.3 and fi A = 0.7 (e.g. Perlmutter et al. 1999). This 
cosmology produces intermediate results to the two cosmolo- 
gies considered here. 



6.1 Implication of the adopted template 

It is worth investigating what happens to the 850-/im flux 
density of M82 as it is placed at progressively higher red- 
shifts, assuming the chosen dust template of = 1.5 and 
Tdust = 40 K. 

Figure ^| demonstrates that in an f2 = 1.0 Universe, 
the 850-/im flux density is almost flat between z = 1 and 
z = 4; the 850-/mi flux density should effectively trace the 
850-^im luminosity. 

In an fl =0.1 Universe, on the other hand, SCUBA 
is less sensitive to z = 4 galaxies than it is to those at 
z = 1. Thus, the higher redshift objects in the sample should 
be more luminous than those at lower redshifts if they are 
detected at the same level. 

The act of increasing (5 and Tdust will effectively make 
SCUBA more sensitive to submillimetre emission at z = 4 
than z — 1. For a given set of fluxes, the luminosities at z — 4 
will be reduced relative to those at z = 1 if /3 and Td us t are 
increased. Likewise, decreasing /3 and Tdust increases the 
inferred 850- luminosities at z = 4 relative to those at 
z = l. 



6.2 hssoiim for the radio galaxy sample 

Table [] gives the 850-/xm luminosities and dust masses for 
the radio galaxy sample. The optically thin, isothermal tem- 
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Table 5. Radio sizes and luminosities for the sample. Column 3 gives the largest angular size (LAS) of the radio source, measured in 
arcseconds. The corresponding linear size (Du near ) and 151-MHz radio luminosity (PisimHz) have been calculated for both fio = 1.0 
and Q = 0.1 (H = 50 kms _1 Mpc -1 ). The radio luminosities were estimated using the radio spectral energy distributions presented 
and referenced in Appendix |a|. The references for the LAS values are given in column 8. Where necessary, the LAS was measured off 
published radio maps. For the majority of the 3CRR sources, the LAS was taken from Blundell et al. (2000), which summarises the radio 
properties of these sources. 



Q = i.o n = o.i 



Source 


Z 


T A Q 
LiAo 


^linear 


1 „ „ T) 

log 1-151MHZ 


^linear 


log -TlSlMHz 


JrtCTS. 






{ ) 


(kpcj 


( \ v Hz sr ) 


(kpcj 


^VVxlz sr ) 




3C277.2 


0.766 


58.0 


473 


27.61 


559 


27.76 


BRR00 


3C340 


0.7754 


46.7 


382 


27.46 


452 


27.60 


BRR00 


3C265 


0.8108 


78.0 


643 


27.86 


768 


28.01 


BRR00 


3C217 


0.8975 


12.0 


101 


27.65 


122 


27.82 


BRR00 


3C356 


1.079 


75.0 


643 


27.92 


809 


28.12 


BRR00 


3C368 


1.132 


7.9 


68 


28.14 


86 


28.35 


BRR00 


3C267 


1.144 


38.0 


327 


28.05 


416 


28.26 


BRR00 


3C324 


1.2063 


10.0 


86 


28.12 


111 


28.34 


BRR00 


3C266 


1.272 


4.5 


39 


28.06 


51 


28.29 


BRR00 


53W069 


1.432 


< 5.1 


< 44 


24.54 


< 59 


24.80 


WvHK84 


4C13.66 


1.45 


6.0 


51 


28.11 


69 


28.36 


BRR00 


3C437 


1.48 


34.4 


294 


28.23 


398 


28.49 


BRR00 


3C241 


1.617 


0.9 


8 


28.30 


11 


28.58 


BRR00 


6C0919+38 


1.65 


10.4 


88 


27.55 


122 


27.84 


NAR92 


3C470 


1.653 


24.0 


203 


28.16 


283 


28.45 


BRR00 


3C322 


1.681 


33.0 


279 


28.19 


390 


28.48 


BRR00 


6C1204+37 


1.779 


51.7 


433 


27.75 


615 


28.06 


LLA95 


3C239 


1.781 


11.2 


94 


28.50 


133 


28.81 


BRR00 


3C294 


1.786 


15.0 


125 


28.39 


179 


28.70 


BRR00 


6C0820+36 


1.86 


23.0 


191 


27.64 


275 


27.95 


LLA95 


6C0905+39 


1.882 


111.0 


920 


27.82 


1330 


28.14 


LGE95 


6C0901+35 


1.904 


2.7 


22 


27.57 


32 


27.90 


NAR92 


5C7.269 


2.218 


7.7 


62 


27.28 


94 


27.64 


BRR00 


4C40.36 


2.265 


4.0 


32 


28.33 


19 


28.69 


CRv097 


MG1744+18 


2.28 


7.5 


60 


28.36 


91 


28.74 


CRv097 


MG1248+11 


2.322 


< 1.2 


< 9 


27.80 


< 15 


28.18 


LBH86 


4C48.48 


2.343 


14.0 


110 


27.95 


171 


28.33 


CRv097 


53W002 


2.39 


0.7 


5 


27.33 


9 


27.71 


WBM91 


6C0930+38 


2.395 


3.7 


29 


27.72 


45 


28.11 


NAR92 


6C1113+34 


2.406 


16.5 


129 


27.68 


201 


28.07 


LLA95 


MG2305+03 


2.457 


3.0 


23 


28.16 


37 


28.55 


SDS99 


3C257 


2.474 


12.0 


93 


28.44 


147 


28.83 


vBS98 


4C23.56 


2.483 


53.0 


411 


28.41 


647 


28.80 


CRv097 


8C1039+68 


2.53 


15.0 


116 


28.00 


183 


28.40 


Lthesis 


MG1016+058 


2.765 


1.3 


10 


27.80 


16 


28.23 


DSD95 


4C24.28 


2.879 


2.3 


17 


28.24 


28 


28.68 


CRv097 


4C28.58 


2.891 


16.2 


119 


28.71 


198 


29.15 


CMvB96 


6C1232+39 


3.221 


7.7 


54 


28.35 


94 


28.82 


NAR92 


6C1159+36 


3.2 


1.2 


8 


28.01 


15 


28.49 


LLA95 


TX1243+036 


3.57 


8.8 


60 


28.49 


107 


29.00 


vOR96 


MG2141+192 


3.592 


8.9 


60 


28.75 


108 


29.26 


CRv097 


6C0032+412 


3.66 


2.3 


15 


28.17 


28 


28.69 


BRE98 


4C60.07 


3.788 


9.0 


59 


28.61 


109 


29.14 


CRv097 


4C41.17 


3.792 


13.1 


86 


28.78 


159 


29.31 


COH94 


8C1435+635 


4.25 


3.9 


2d 


28.42 


47 


28.99 


LMR94 


6C0140+326 


4.41 


2.5 


15 


28.20 


30 


28.78 


RLB96, BRE98 



plate spectrum with /? = 1.5 and Td us t = 40K has been as- 
sumed. The errors on the values of Lsso^m and Md us t have 
been calculated using the errors in the flux density mea- 
surements only; they do not account for the uncertainties in 
cosmology, /3, or T dus t. 

The dust mass of each galaxy in the sample was calcu- 
lated from the 850-pm observation using: 



Mdust = -r— r j= r (1) 

(1 + ZJ Ripest tSu TSBt \lduat) 

where S b s is the observed flux-density, u res t is the cor- 
responding rest-frame frequency, Dl is the luminosity dis- 
tance to the source, K Vrcst is the mass absorption coeffi- 
cient of the dust at iVest, and B Vrl!Bt (Tdust) is the intensity 
of a blackbody at v res t assuming isothermal emission from 
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dust grains at a temperature Tdust- The mass absorption 
coefficient k is poorly co nstrained, with publ ished values of 
ft(800^m): 0.04 m 2 kg _1 ( [Draine fc Lee 1984| ), 12 nAg' 1 
(Chini, Kriigel, & Kreysa T986), 0.15 m^kg^ 1 flHildebrand 
1983|)71ind 0.3 m 2 kg _1 flVlathis fc Whiffen 1989[ ). These 
estimates c an be extrapolated to other wavelengths using 
k oc A _/3 (Chini et al. 1986). The intermediate value of 
Ac(800^im)=0.15 m kg~ will be assumed here. A different 
choice will not affect the relative dust masses of the sample, 
and can only alter the absolute dust masses by a factor of 
~ 4 at most ( Hughes et al. 1997 ) . 

In Figure [ij the 850- /im flux density is plotted against 
redshift for the radio galaxy sample. Plots of Lsso^m vs. 
z in Q.o = 1.0 and Q. a = 0.1 Universes are also shown. As 
expected from Figure^, in the tt a = 1.0 case the flux density 
traces the luminosity between z = 1 and z = 4 , and the two 
plots look the same. For fi D = 0.1, the luminosity of the 
z ~ 4 galaxies is higher relative to those at z ~ 1 , as 
SCUBA is less sensitive to the high-redshift objects. 

The small dynamic range in the measured values of 
Lssofjm for this study are limited by the sensitivity of 
SCUBA and the characteristic dust mass in high-redshift 
objects. Note the detections of L^so^™ appear at roughly the 
same level as the Lgso^m upper limits; most of the galaxies 
in the sample have been observed to the same depth in 850- 
/im luminosity. Coupled with the fact that the upper limits 
are clustered at the low-redshift end of the sample, this in- 
dicates that the high-redshift radio galaxies are intrinsically 
brighter in the submillimetre than the lower-redshift radio 
galaxies. This evidence for evolution of the radio galaxy pop- 
ulation with redshift will now be analysed in detail. 



7 EVIDENCE FOR EVOLUTION 

7.1 Survival analysis: correlations between radio 
properties, submillimetre properties, and 
redshift 

A logical approach to finding evidence which either sup- 
ports or rejects the evolution scenario is to run correlation 
tests between the radio properties, submillimetre properties, 
and redshifts of the sample. For example, a correlation be- 
tween L850fim and redshift, in the absence of a correlation 
between Lsso^m and Pisimhz, will imply that radio galaxies 
undergo cosmological evolution of their submillimetre emis- 
sion. Conversely, a significant correlation between Lsso/im 
and Pisimhz in the absence of a correlation between Lg5o Mm 
and redshift, would indicate that the submillimetre proper- 
ties of radio galaxies are more closely tied to the radio prop- 
erties of the galaxy (perhaps, for example, via host mass) 
than to cosmological epoch. 

The obvious correlations to consider are between 
Lssoum, Pisimhz, and z. It is also instructive to investi- 
gate whether Pisimhz correlates with Du near (the linear size 
of the radio source), and whether the radio spectral index 
aradio correlates with either L 850Al m, Pisimhz, or z. 

As mentioned previously, for radio sources selected from 
flux-limited samples, the 151-MHz radio luminosity should 
be a good indicator of the intrinsic power Q of the radio jet. 
The size of a radio source depends on its age, the medium 
it is expanding into, and the jet power. Again, consider- 
ing sources from flux-limited samples, Willott et al. ( 199E ) 
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Figure 4. Scatter plots of 850-fim flux density and (rest-frame) 
luminosity against redshift for the radio galaxy sample. Lsso^m 
has been calculated assuming j3 = 1.5 and T(j us t = 40K, and is 
shown for both f2 = 1.0 and Q a = 0.1 (H is assumed to be 
50 kms-iMpc" 1 ). 



showed that the dependence of Dunear on jet power is very 
weak; Dn ne ar correlates most strongly with the age of the 
radio source. Although radio luminosities are expected to 
decline throughout the lifetime of radio sources (Kaiser, 
Dennett-Thorpe & Alexander 1997; Blundell et al. 1999; 
Blundell & Rawlings 1999), selection effects conspire to en- 
sure a very weak anti-correlation of radio luminosity with 
D(i„ear for sourc es selected from a nu mber of separate flux- 
limited samples (Blundell et al. 1999). Although our sample 



of SCUBA targets is not purely a combination of flux- limited 
samples, it is dominated by 3C/3CRR and 6CE targets so 
no strong Pisimhz-D links are expected. 

The criterion that a ra di is large is often used to identify 
high redshift sources in large radio surveys. This selection is 
based on the empirical correlation observed between a ra dio 
and radio luminosity (e.g. Veron, Veron & Witzel 1972). In 
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Table 6. Rest-frame 850-/xm luminosities (Lg5o Mm ) and dust masses (Md us t) for the radio galaxy sample. Optically thin, isothermal 
greybody emission has been assumed, with (3 = 1.5 and T duat = 40K. Errors have been estimated using only the standard errors on the 
flux density measurements; 2<r upper limits are given for undetected (S/N < 2) sources. Lsso^m and M dust have been calculated for 
three different cosmologies: Sl o =1.0, H o =50 kms - 1 Mpc _1 ; Q o =0.1, H o =50 kms - ^Mpc -1 ; f2 o =0.1, H =67 kms -1 Mpc — 1 . Note, column 
3 indicates whether the source belongs to either or the subsets defined in Section [7.1. 



Q o =1.0, H o =50 C o =0.1, H o =50 n o =0.1, H =67 

Source z Subset log L 8 50^m log M dust log L 850Mm log M dust log L g50Mm log M duat 

(WHz-isr" 1 ) (M s ) (WHz-isr" 1 ) (M Q ) (WHz^sr 1 ) (M Q ) 



3C277.2 


0.766 




< 


22.70 


< 


8.18 


< 22.85 


< 8.32 


< 22.59 


< 


8.07 


3C340 


0.775 




< 


22.63 


< 


8.10 


< 22.77 


< 8.25 


< 22.52 


< 


8.00 


3C265 


0.811 




< 


22.69 


< 


8.16 


< 22.84 


< 8.32 


< 22.59 


< 


8.06 


3C217 


0.897 




< 


22.63 


< 


8.10 


< 22.80 


< 8.27 


< 22.54 


< 


8.02 


3C356 


1.079 




< 


22.84 


< 


8.31 


< 23.04 


< 8.51 


< 22.78 


< 


8.26 


3C368 


1.132 






23.00 




8.47 


23.20 


8.68 


22.95 




8.43 


3C267 


1.144 




< 


22.71 


< 


8.19 


< 22.92 


< 8.40 


< 22.67 


< 


8.14 


3C324 


1.206 




< 


22.84 


< 


8.32 


< 23.06 


< 8.54 


< 22.81 


< 


8.28 


3C266 


1.272 




< 


22.85 


< 


8.32 


< 23.08 


< 8.55 


< 22.82 


< 


8.30 


53W069 


1.432 




< 


22.75 


< 


8.22 


< 23.00 


< 8.48 


< 22.75 


< 


8.22 


4C13.66 


1.45 


A,B 




22.85 




8.32 


23.11 


8.58 


22.85 




8.33 


3C437 


1.48 


A,B 


< 


22.72 


< 


8.20 


< 22.98 


< 8.46 


< 22.73 


< 


8.20 


3C241 


1.617 


A,B 


< 


22.79 


< 


8.26 


< 23.07 


< 8.55 


< 22.82 


< 


8.29 


6C0919+38 


1.65 




< 


22.75 


< 


8.22 


< 23.03 


< 8.51 


< 22.78 


< 


8.25 


3C470 


1.653 


A,B 


< 


23.21 


< 


8.69 


< 23.50 


< 8.98 


< 23.25 


< 


8.72 


3C322 


1.681 


A,B 


< 


22.75 


< 


8.22 


< 23.04 


< 8.52 


< 22.79 


< 


8.26 


6C1204+37 


1.779 




< 


22.81 


< 


8.29 


< 23.12 


< 8.60 


< 22.87 


< 


8.34 


3C239 


1.781 


A,B 


< 


22.72 


< 


8.19 


< 23.02 


< 8.50 


< 22.77 


< 


8.25 


3C294 


1.786 


A,B 


< 


22.61 


< 


8.09 


< 22.92 


< 8.39 


< 22.66 


< 


8.14 


6C0820+36 


1.86 




< 


22.99 


< 


8.46 


< 23.30 


< 8.78 


< 23.05 


< 


8.53 


6C0905+39 


1.882 






22.96 




8.43 


23.28 


8.75 


23.02 




8.50 


6C0901+35 


1.904 




< 


22.77 


< 


8.25 


< 23.10 


< 8.57 


< 22.84 


< 


8.32 


5C7.269 


2.218 




< 


22.93 


< 


8.40 


< 23.29 


< 8.77 


< 23.04 


< 


8.51 


4C40.36 


2.265 


A,B 
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22.83 


< 


8.31 


< 23.20 
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< 


8.42 


MG1744+18 


2.28 


A,B 


< 


22.70 


< 


8.18 


< 23.08 


< 8.55 


< 22.82 


< 


8.30 


MG1248+11 


2.322 




< 


22.88 


< 8.36 


< 23.26 


< 8.73 


< 23.00 


< 


8.48 


4C48.48 


2.343 


A,B 




23.06 




8.54 


23.44 


8.92 


23.19 




8.66 


53W002 


2.39 




< 


22.89 


< 


8.36 


< 23.27 


< 8.75 


< 23.02 


< 


8.49 


6C0930+38 


2.395 




< 


22.69 


< 


8.16 


< 23.07 


< 8.55 


< 22.82 


< 


8.29 


6C1113+34 


2.406 




< 


22.74 


< 


8.22 


< 23.13 


< 8.60 


< 22.87 


< 


8.35 


MG2305+03 


2.457 


A,B 


< 


22.71 


< 


8.18 


< 23.10 


< 8.57 


< 22.84 


< 


8.32 


3C257 


2.474 


A,B 


< 


23.17 


< 


8.64 


< 23.56 


< 9.04 


< 23.31 


< 


8.78 


4C23.56 


2.483 


A,B 


< 


22.94 


< 


8.42 


< 23.34 


< 8.81 


< 23.08 


< 


8.56 


8C1039+68 


2.53 


A,B 


< 


22.68 


< 


8.15 


< 23.08 


< 8.55 


< 22.82 


< 


8.30 


MG1016+058 


2.765 


B 




22.73 




8.20 


23.16 


8.63 


22.90 




8.38 


4C24.28 


2.879 


A,B 




22.73 




8.21 


23.17 


8.65 


22.92 




8.39 


4C28.58 


2.891 






22.95 




8.42 


23.39 


8.86 


23.13 




8.61 


6C1232+39 


3.221 


A,B 




22.93 




8.40 


23.40 


8.88 


23.15 




8.62 


6C1159+36 


3.2 


A,B 


< 


22.83 


< 


8.31 


< 23.31 


< 8.78 


< 23.06 


< 


8.53 


TX1243+036 


3.57 


A 


< 


22.97 


< 


8.45 


< 23.48 


< 8.95 


< 23.22 


< 


8.70 


MG2141+192 


3.592 






22.99 




8.47 


23.50 


8.98 


23.25 




8.72 


6C0032+412 


3.66 


A,B 




22.71 




8.18 


23.22 


8.70 


22.97 




8.45 


4C60.07 


3.788 






23.56 




9.03 


24.08 


9.56 


23.83 




9.31 


4C41.17 


3.792 






23.40 




8.88 


23.93 


9.41 


23.68 




9.16 


8C1435+635 


4.25 


A 




23.20 




8.68 


23.77 


9.25 


23.52 




8.99 


6C0140+326 


4.41 


A,B 




22.83 




8.30 


23.41 


8.89 


23.16 




8.63 



a flux-limited survey, high values of a ra dio will often corre- 
spond to a large intrinsic luminosity and hence to a high 
redshift. 

However, it is possible to envisage models in which 
Q-radio is related closely to star-formation activity (e.g. Lilly 
1989). This is relevant to the work presented here, as sev- 
eral of the high-redshift galaxies in the sample were origi- 
nally identified in these ultra-steep spectrum searches (e.g. 



Chambers, Miley & van Breugel 1988, 1990; Chambers et al. 
1996; Stern et al. 1999). Thus, if Lsso^m appears correlated 
with redshift, it could be a manifestation of how the sources 
were selected in the first place, a possibility which can be 
tested by investigating whether a correlation exists between 
Lsso/im and a ra dio- Note, for the correlation analysis, a ra dio 
was calculated at 151 MHz. 
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Table 7. Results of the survival analysis correlation tests applied to the entire sample. The significance (P) is the probability of the null 
hypothesis (no correlation) being true. If all three tests yield P < 5 per cent, the variables are taken to be correlated. If only some of 
the tests yield P < 5 per cent, evidence for a correlation is taken to be uncertain. 



Cosmology 


Variable 


Percentage of 




Significance (P) 


Correlation 




Dependent 


Independent 


Data Censored 


Uox 


rvcnaaii 


Spearman 


Present? 


- 


SsSO^tm 


z 


70% 


0.00% 


0.02% 


0.01% 


YES 




n, . 

Linear 


l°g(Pl51Mifz) 


4% 


22.48% 


CO RQ% 

oz.oy /o 


(1.1/ /o 


NO 




io g(r 151MHz ) 


z 


u/o 


U.U1/0 


U.Uo/o 


U.U f /o 


VTTQ 


n a = 1.0 




io gt^l51MJfzJ 


70% 


2.65% 


U.zl /o 


U.D4/0 


YES 


log(L 85 o Mm ) 


Z 


70% 


0.00% 


0.19% 


0.08% 


YES 




log(Lg50/jm) 




70% 


0.31% 


0.36% 


1.73% 


YES 




^radio 


log(Pi5iMJJ 2 ) 


0% 


0.64% 


1.00% 


1.36% 


YES 




^radio 


Z 


0% 


4.18% 


13.94% 


10.50% 


MAYBE 




^linear 


log(PltlA/JJ 2 ) 


4% 


32.29% 


67.59% 


50.22% 


NO 




l°g(Pl51MHz) 


Z 


0% 


0.00% 


0.00% 


0.00% 


YES 


n = o.i 


log(L850^m) 


iog(Pi5iMH^) 


70% 


0.53% 


0.01% 


0.08% 


YES 


log(L g5 o Mm ) 


z 


70% 


0.00% 


0.00% 


0.00% 


YES 




log(L 850Mm ) 


C^radio 


70% 


0.18% 


0.67% 


1.88% 


YES 






log(Pi5iMJJ 2 ) 


0% 


0.44% 


1.61% 


1.47% 


YES 




^radio 


z 


0% 


4.18% 


13.94% 


10.50% 


MAYBE 


Total Number of Galaxies: 46 



7.1.1 Survival analysis 

The submillimetre dataset for the radio galaxy sample con- 
tains a large number of upper limits, or censored data. Run- 
ning correlation tests on the dataset is not a trivial matter, 
as standard statistical theory is not equipped to handle this 
kind of information. 

Statisticians have extensively studied the problem of 
censored datasets, and a means for successfully handling 
them has been known for decades. It is called survival anal- 
ysis and it allows censored data samples to be treated in a 
meaningful manner with minimal loss of information - the 
values of the upper /lower limits are formally taken into ac- 
count when calculating correlation coefficients and other test 
statistics. 

The application of survival analysis to astronomical 
data is outlined in two papers: Feigelson & Nelson 1985 (Pa- 
per I) and Isobe, Feigelson, & Nelson 1986 (Paper II). Paper 
I deals with univariate problems. Paper II considers bivariate 
techniques to characterise the correlation and line-regression 
between two variables. The test statistics described in these 
papers are available in the software package ASURV Rev. 1.1 
(Isobe & Feigelson 1990; LaValley, Isobe & Feigelson 1992). 

ASURV contains three different correlation tests - Cox's 
Proportional Hazard Model, Generalized Kendall's Tau, and 
Generalized Spearman's Rho. All three methods test the null 
hypothesis that no correlation is present. The ASURV rou- 
tines compute the probability of the null hypothesis being 
true, P, also referred to as the significance of the correla- 
tion. The convention to be adopted here is that X and Y are 
correlated if P < 5 per cent. A correlation will be treated as 
firm if all three tests yield P < 5 per cent; the existence of 
a correlation will be treated as possible, but unconfirmed, if 
only some of the tests indicate P < 5 per cent. 

Note, the Cox test only allows censoring in the depen- 
dent variable; the Kendall and Spearman tests allow censor- 



ing in both the dependent and the independent variables. 
The Spearman technique has not been properly tested yet. 
It is known to give misleading results if < 30. 



7.1.2 Correlations for the entire dataset 

Initially, the entire sample was used to run the correlation 
tests. The results of this analysis are given in Table 

Figure ^ predicts that the 850-/im flux density will mir- 
ror the 850-/im luminosity if Q. a = l.Q. This is reflected by 
the S850Mm-z and Lsso/xm-z correlations having similar sig- 
nificances. As expected, decreasing Q to 0.1 increases the 
strength of the Lgso/jm-z correlation. 

The absence of a correlation between Pisimhz and 
Ounear is indicative of a c ombin ed sample in which, as dis- 
cussed by Willott et al. (1999), Pisimhz can be taken to 
crudely track jet power, and D;;„ ear to crudely track source 
age. 

It is troublesome, but not unexpected, that Lsso^m, 
Pisimhz, and z are all correlated with each other. The 
strength of these correlations is greater in a low-f2 Uni- 
verse. The strong correlation between Pisimhz and z could 
be responsible for either the Lsso/jm-z correlation if Lsso/jm 
actually correlates with Pisimhz, or for the L 8 50Mm-PisiMHz 
correlation if Lg5o Mm actually correlates with z. Thus, if the 
effects of cosmological evolution and radio power are to be 
disentangled, the entire dataset cannot be used. 

L85ivm also appears correlated with a ra dio. However, 
this may be a remnant of the fact that a ra dio and Lsso/jm 
are both correlated with Pisimhz and z. Again, if the entire 
dataset is used, discriminating between different effects is 
difficult. 

It was decided to choose a strip in the P-z plane that 
covered a large range of redshifts in a single narrow band of 
radio power. In the absence of a correlation between Pisimhz 
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Table 8. Results of the survival analysis correlation tests applied to Subset A. The significance (P) is the probability of the null 
hypothesis (no correlation) being true. The generalized Spearman test has not been applied as it breaks down if the sample size is less 
than 30. If both Cox and Kendall yield P < 5 per cent, the variables are taken to be correlated. If only one of them yields P < 5 per 
cent, evidence for a correlation is taken to be uncertain. 



Cosmology 


Variable 


Percentage of 


Si! 


^nificance (P) 


Correlation 




Dependent 


Independent 


Data Censored 


Cox 


Kendall Spearman 


Present? 




S850m"i 


z 


67% 


0.29% 


4.05% 


YES 




^linear 


log(Pi5iA/tf s ) 


0% 


63.20% 


34.86% 


NO 




log(Pi5iAIff z ) 


Z 


0% 


74.43% 


50.65% 


NO 


n = 1.0 


log(L g5 o Mm ) 


log(Pi5iA/tfz) 


67% 


50.52% 


77.71% 


NO 


log(L g5 o^ m ) 


Z 


67% 


1.91% 


15.70% 


MAYBE 




log(L 850Mm ) 


C^radio 


67% 


70.37% 


85.01% 


NO 




CX-radio 


log(Pl51AfH Z ) 


0% 


21.03% 


16.40% 


NO 




^radio 


z 


0% 


21.22% 


18.32% 


NO 



Total Number of Galaxies: 21 



Table 9. Results of the survival analysis correlation tests applied to Subset B. The significance (P) is the probability of the null 
hypothesis (no correlation) being true. The generalized Spearman test has not been applied as it breaks down if the sample size is less 
than 30. If both Cox and Kendall yield P < 5 per cent, the variables are taken to be correlated. If only one of them yields P < 5 per 
cent, evidence for a correlation is taken to be uncertain. 



Cosmology Variable 

Dependent Independent 


Percentage of 
Data Censored 


Si! 

Cox 


jnificance (P) 
Kendall Spearman 


Correlation 
Present? 


S850Mm 


z 


65% 


0.81% 


5.87% 


MAYBE 


^linear 


l°g(Pl51A/ff z ) 


0% 


84.35% 


49.54% 


NO 


log(Pi5ijv/Hz) 
log(Lg50 Mm ) 

a c = o.i io g (L 850Mm ) 


z 

log(Pi5i A / ffz ) 

z 


0% 
65% 
65% 


34.86% 
79.95% 
0.13% 


27.00% 
96.22% 
2.01% 


NO 
NO 
YES 


log(L850,j m ) 


O^radio 


65% 


64.95% 


92.43% 


NO 



a radio log(Pi5iA/ Hz ) 0% 0.39% 2.49% YES 

a radio z 0% 20.64% 41.66% NO 

Total Number of Galaxies: 20 



and z, it could then be determined whether Lsso^m truly 
correlates with z or Pisimhz- The strip was chosen to satisfy 
the following criteria: 

(i) To cover the largest range in redshift possible whilst 
still maintaining even coverage of radio power. This auto- 
matically imposed the restriction that z > 1.4 due to the 3C 
survey having a lack of low-redshift, very high radio-power 
sources. 

(ii) To contain at least 20 objects to satisfy the require- 
ments of the survival analysis techniques. 

Initially a strip was chosen using the P-z plane for 
Qo~l-0, hereafter referred to as Subset A. A different but 
similar strip was chosen from the P-z plane for £7 D =0. 1 . It 
will be referred to as Subset B. The sources included in each 
subset are listed in Table || 

It is important to note that there is nothing special 
about the strips. They have merely been chosen to help dis- 
entangle the various influences affecting Lssoum • The results 
found for the subsets should apply to the entire sample and 
to radio galaxies in general. 



7.1.3 Correlations for Subset A 

Subset A, shown in Figure ^, contains 21 radio galaxies and 
is defined by z > 1.4 and 27.9 < log(Pi 5 i M Hz) < 28.5 for 
fi o = 1.0, Ho = 50 kms^Mpc- 1 . 

Table ^| contains the results of the survival analysis ap- 
plied to Subset A, and Figure |E| contains the corresponding 
scatter plots. 

For f2 o =1.0, the existence of a correlation between 
Lsso,™ and redshift is hinted at, but is not confirmed by 
all of the correlation tests. The important point, however, 
is that L850fim definitely fails to correlate with Pisimhz and 
otradio- The probability of 'no correlation' between Lsso/jm 
and z is much lower than for Lgso/jm-PisiMHz or LssOfjm- 

Survival analysis could also be applied to Subset A for 
the f2 o =0.1 case. However, for Subset A, Pisimhz and z are 
correlated if f2 o =0.1. This is precisely what was trying to 
be avoided by selecting the subset in the first place. Thus, 
a separate subsample will need to be chosen to study the 
statistics for a low-S7 Universe. 
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Figure 5. Scatter plots corresponding to the correlations investigated for Subset A. 



7.1.4 Correlations for Subset B 

Subset B, shown in Figure ^, contains 20 radio galaxies and 
is denned by z > 1.4 and 28.2 < log(Pi 5 i M Hz) < 28.9 for 
Qa=0.1, Ho = 50 kms^Mpc" 1 . 

Table ^| contains the results of the survival analysis ap- 
plied to Subset B, and Figure [j] contains the corresponding 
scatter plots. 

The results of the survival analysis are more clear cut 



in a low-f2 Universe. The suspected correlation between 
Lsso,™ and z is confirmed as significant by all of the corre- 
lation tests. Neither z nor Lsso^m correlates with Pisimhz- 
This indicates that the important correlation for Lssopm is 
most likely to be with redshift, not radio luminosity. 

For Subset B, Lsso/jm does not correlate with a radio or 
Pisimhz. However, a radio and Pisimhz correlate with each 
other. This implies that the correlation between a ra dio and 
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Figure 6. Scatter plots corresponding to the correlations investigated for Subset B. 



LssoMm seen for the entire dataset is an artifact of their mu- 
tual correlation with Pisimhz- 

7.1.5 Considering the dust emission template 

The analysis described so far assumes a template for dust 
emission with (3 = 1.5 and Td us t = 40K. However, what if 
this template was incorrect? Would the results of the analy- 
sis be affected if the true values of /3 and Tdust were different? 



It has already been observed that increasing f3 and Tdust 
makes SCUBA more sensitive to submillimetre emission at 
z = 4 than it is at z = 1. Given the measured fluxes, the 
850-pim luminosities at z — 3 — 4 will be reduced relative 
to those at z = 1 — 2 if j3 and Tdust are increased. As a 
result, correlations between Lgso Mm and z should be weak- 
ened. Likewise, decreasing f3 and Td us t should strengthen 
Lgso/jm-z correlations. 

Changing f3 and Td us t will not alter Lssopm in a radio- 
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Table 10. Summary of the Lssopm-z correlation analysis for the 
two cosmologies considered and for reasonable values of j3. / in- 
dicates the survival analysis found a correlation between Lssojjm 
and z, X indicates the survival analysis found no such correlation, 
and ? indicates that some, but not all, of the correlation tests 
found a correlation to be present. 



a, = 1.0 Via = 0.1 







Entire 
Datasct 


Subset 
A 


Entire 
Datasct 


Subset 
B 


= 


1.0 


/ 


/ 


/ 


/ 


(3 = 


1.5 


/ 


? 


/ 


/ 


& = 


2.0 


/ 




/ 


/ 



Table 11. Detection rate of radio galaxies at 850 fim for z < 2.5 
and z > 2.5. For each dataset and redshift interval, the number 
of galaxies present is given in parenthesis. Note, 6C1159+36 may 
be at a lower redshift than previously published, which would 
increase the detection rate at z > 2.5. 





Detection Rate 


Dataset 


z < 2.5 


z > 2.5 


Subset A 


15% (13) 


63% (8) 


Subset B 


15% (13) 


71% (7) 


All Data 


12% (33) 


77% (13) 



power or radio-spectral index dependent manner. If the dust 
template is wrong, the correlations between Lsso^m and 
P i5iMHz /a ra dio should not be affected in a systematic way. 

In order to investigate this fully, the survival analysis 
was re-run using several different combinations of /3 and 
Tdust- To examine the effect of changing f3, Td us t was held at 
40 K while /3 was increased to 2.0, and decreased to 1.0. Like- 
wise, P was held at 1.5 while Tdust was changed to 100 K, 
70 K, and 20 K. 

The observed effects on the results of the correlation 
tests can be summarised as follows: 

(i) L850mwi-z: Increasing [3 or Tdust destroys the marginal 
correlation between Lgso Mm and z for Subset A. For the en- 
tire dataset (both cosmologies) and for Subset B, the cor- 
relation is evident regardless of the values of /3 and Tdust- 
Decreasing j3 and Td us t confirms the possible Lgso^m-z cor- 
relation for Subset A. 

(ii) L 8 50 M m-Pi5iMHz: Changing (3 and T dus t does not de- 
stroy the Lg5o M m-Pi5iMHz correlation for the entire dataset, 
nor does it create a correlation between Lsso^m and Pisimhz 
for Subsets A or B. 

(iii) L 85 o^ m -a radlo : Changing (3 and T dus t does not de- 
stroy the ljS50fj,m-a ra dio correlation for the entire dataset, 
nor does it create a correlation between Lsso^m and a ra dio 
for Subsets A or B. 



7.1.6 Summary of correlation analysis 

For Subsets A and B, survival analysis finds Lgsoum to be 
correlated with redshift unless the highest values of /3 and Q. a 
are chosen (Table |l(J). Furthermore, Lsso^m is not found to 
correlate with either Pisimhz or a ra di for any of the tested 
cosmologies or dust emission templates. Even for high values 
of /3 and Q , there is evidence that Lgso^m is dependent on 
redshift instead of on the radio properties of the source: 
the probability of 'no correlation' is substantially smaller 
for Lgsoum and z than for Lgso^m and Pisimhz or a radio . 

This all suggests that the correlation between Lsso^m 
and redshift for the entire dataset is real, and is not an ar- 
tifact of their mutual correlation to Pisimhz (although this 
may affect the strength of the observed Lg5o Mm -2 correla- 
tion). 

Given the low detection threshold adopted in this paper, 
it is worth re-iterating that a detection threshold of S/N> 3 
(with upper limits taken at the 3<r level) yields the same 
result. 



7.2 Further investigations of the relationship 
between Lssonm. and redshift 

SCUBA has had more success detecting the higher-redshift 
galaxies in the sample. This is shown in Table jn], which 
display the 850-pm detection rate for z < 2.5 and z > 2.5. 
At 850 fj,m, SCUBA is almost equally sensitive to all of the 
galaxies in the sample. The striking result that almost none 
of the galaxies are detected if z < 2.5, whereas almost all of 
them are detected if z > 2.5, suggests a trend of increasing 
submillimetre luminosity with redshift. Note, the detection 
rates observed for Subset A and Subset B closely mirror 
those found for the entire sample, providing further evidence 
for the cosmological evolution of the radio galaxy population 
irrespective of the radio source. 

There is an alternative to survival analysis for investi- 
gating the relationship between Lsso^m and z. Instead of 
running a correlation test which only takes the strength 
of the detection or upper-limit into account, the 850-/im 
luminosities can be binned in redshift to calculate the 
typical submillimetre luminosity as a function of redshift. 
This method considers both the signal and error for each 
source and is independent of the precise choice of detec- 
tion threshold. Note, for the following discussion, H has 
been taken to be 67 kms _1 Mpc _1 for fl a = 0.1; this should 
give a comparable measure of luminosity to the fi = 1.0, 
Ho = 50 kms~ 1 Mpc _1 cosmology. 

Four redshift bins have been chosen which are ade- 
quately sampled by the dataset and are centered as near as 
possible on z = 1, 2, 3, 4: 0.45 < z < 1.45, 1.45 < z < 2.45, 
2.45 < z < 3.45, and 3.45 < z < 4.45. The weighted-mean 



( Bevington_ 1969 ) 850-/mi luminosity, <L85o^m>, is plotted 
in Figure ^ for the entire dataset and Subset A assuming 
an fi = 1.0 cosmology and for the entire dataset and Sub- 
set B assuming an fi =0.1 cosmology. Curves of the form 
L = a(l + z) b have been fit to <L850/im> in order to help 
quantify the strength of the 850-/im luminosity evolution. 

There are several important points to be made from 
Figure |?J 

(i) For the entire dataset, there are clear steps in 
<L850/jm> between each redshift bin. If this trend is char- 
acterised by the function L = a(l + z) b , nonlinear least- 
squares fitting indicates that b = 4.58 ± 0.52 if S7 o = 1.0, and 
b = 5.52 ±0.53 for Q o =0.1. 

(ii) For Subset A, there are clear steps in <L850/jm> be- 
tween the z = 2 and z = 3 bins, and between the 2 = 3 
and z = 4 bins. If this trend is characterised by the function 
L = a(l + z) b , nonlinear least-squares fitting indicates that 
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Figure 7. Weighted-mean 850-^tm luminosity, <Lg5o^m>, binned in redshift. <Lgso (im > has been plotted at the midpoint of each bin, 
which at most varies by 0.2 from the average redshift of all the objects in the bin. The error bars only represent the statistical errors on the 
observations; they do not account for uncertainties in either the cosmology or the dust template. Curves of the form L = a(l + z) b have 
been fit to the data using nonlinear least-squares fitting, and the horizontal line indicates the origin (or Lsso/jm = 0.0). The upper plot 
displays an f2 = 1.0, H = 50 kms - 1 Mpc~ 1 cosmology, and <Lg5o Mm > is plotted for both the entire dataset and Subset A (diamond 
symbols). The solid curve is the best-fit for the entire dataset, with a = 0.010 ± 0.008, b = 4.58 ± 0.52, x 2 = 1-23. The dashed curve is 
the best-fit for Subset A, with a = 0.029 ± 0.034, b = 3.67 ± 0.76, x 2 = 0.29. Similarly, the lower plot displays the entire dataset and 
Subset B (diamond symbols) for an f2 = 0.1, H = 67 kms _1 Mpc _1 cosmology. The solid curve is the best-fit for the entire dataset, 
with a = 0.004 ± 0.003, b = 5.52 ± 0.53, x 2 = 1-58. The dashed curve is the best-fit for Subset B, with a = 0.070 ± 0.099, b = 3.23 ± 1.01, 
X 2 = 0.32. Subsets A and B do not contain any objects in the z = 1 bin. 



b = 3.67 ±0.76. Subset A does not contain any objects in 
the z = 1 bin. 

(iii) For Subset B, there is a clear step in <L850/jm> be- 
tween the z — 2 and z = 3 bins. There may also be an in- 
crease in <Li850/«n> from z = 3 to z — 4, although the error 
bars do not preclude <L85o Mm > flattening off between these 
two redshifts. If the model L = a(l + z) b is assumed, nonlin- 
ear least-squares fitting indicates that b = 3.23 ± 1.01. Sub- 
set B does not contain any objects in the z — 1 bin. 

(iv) For Subsets A and B, the correlation analysis has re- 
vealed that L850(jm does not depend on the radio properties 
of the source. These subsamples are not special, they have 
been chosen to cover an even spread in radio luminosity over 
a range of redshifts; their selection has been dictated only by 
the deficiencies of the original radio surveys. Thus, the evolu- 



tion of <L85o,jm > with redshift observed for these two sub- 
samples, which is consistent with < L >~ (1 + z) 3 , should 
apply to all radio galaxies. This is supported by the fact that 
for the z — 2 and z = 3 bins, the value of <L85o f ™ > for the 
subset mimics the value of <L85o,jm> for the entire dataset. 
At z ~ 4, however, some of the most extreme objects in the 
Universe have been observed (e.g. 4C60.07 and 4C41.17), 
which are simultaneously the most radio-luminous and the 
most submillimetre-luminous, residing at the highest red- 
shifts. These galaxies fall outside the subsets which have 
been constructed to exclude radio-luminosity bias. Thus at 
z ~ 4, <L85o M m> for the entire sample far exceeds that of 
the subsets. 

As for the correlation analysis, it is important to in- 
vestigate the robustness of this result to variations in the 
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dust template. Thus, (3 was changed to 1.0 and 2.0, whilst 
Tdust was fixed at 40K. Likewise, holding f3 constant at 1.5, 
Tdust was changed to 20K, 70K, and 100K. Altering the dust 
template has a larger effect on the higher redshift sources. 
Thus, if /3 and Tdust are increased, the trend of increasing 
<Lg5o,jm> with redshift is weakened; if they are decreased 
the trend is strengthened. More specifically, it was found 
that: 

(i) For the entire dataset, there are always clear steps in 
<Lg5o Mm > between all redshift bins, regardless of the values 
of (3 and T dust . 

(ii) For Subset A, there is always a clear step in 
<Lgso Mm > between z = 2 and z = 3 and between z — 3 
and z — 4, regardless of the values of j3 and Tdust ■ 

(iii) For Subset B, there is always a clear step in 
<Lgso M m> between z = 2 and z = 3. The possible step 
between z — 3 and 2 = 4 becomes clear (i.e. the error bars 
do not allow <Lg5o Mm > to remain flat) if j3 is decreased to 
1.0 or if Tdust is decreased to 20K. 

(iv) For the z = 2 and z = 3 bins, <Lg5o M m> for the en- 
tire dataset always closely mimics <Lg5o M m> for the subsets, 
regardless of the values of /3 and Tdust. 

(v) Adopting the model L = a(l + z) b , the luminosity 
evolution observed For Subsets A and B upon changing j3 
and Tdust is consistent with 'b' lying between 2.5 and 5.0. 

Binning the data has been able to confirm the relation- 
ship between Lgso^m and z that was initially suggested by 
survival analysis. It is therefore possible that binning may 
also reveal a relationship between Lg5o Mm and Pisimhz or 

Qradio ■ 

It is only necessary to check this for Subsets A and 
B, as Lg5o, lm is known to be strongly correlated with both 
Pisimhz and 

Otradio 

for the entire dataset. If Lgso^m remains 
independent of the radio source for Subsets A and B, the 
evidence for evolution with redshift will be strengthened. 

Binning Lgso^m against radio power and radio spectral 
index fails to reveal a dependence of Lgso,™ on the radio 
source. Altering the dust template does not change this re- 
sult. 



7.3 Addressing potential remaining selection 
effects 

7.3.1 Potential contamination by obscured quasars 

Careful scrutiny of the P-z plane defined by combining all 
the known high-redshift radio sources shows that there is 
a lack of high-redshift, high-radio power quasars, although 
this is at least partly due to the selection effects inherent 



in stee p-spectrum-selected radio samples (Jarvis et al., in 
preparation). This is also the area of the P-z plane where the 
850-^m detections of the radio-galaxy sample are clustered. 

According to Unified Schemes, radio-loud quasars and 
radio galaxies are identical phenomena: the former is viewed 
within the opening angle of an obscuring dust torus allowing 
the quasar to be seen, the latter is viewed outside the open- 
ing angle of the dust torus which hides the quasar from view 
(e.g. Antonucci & Barvainis 1990 and references therein). 

It can be argued that at higher redshifts, there is an 
increased chance of dust hiding the optically-bright nucleus 
of a radio-loud quasar, even though the radio-loud quasar is 



viewed within the opening angle of the dust torus. Thus, po- 
tentially, several high-redshift radio-loud quasars may have 
been mis-classified as radio galaxies. Therefore, the worry for 
the present study is that our radio galaxy sample could be- 
come progressively more contaminated by obscured quasars 
with increasing redshift. 

Of course, within unified schemes all powerful radio 
galaxies contain a buried quasar. Thus the point of issue 
for the present study is whether the submillimetre emission 
of a powerful radio source is orientation dependent. If so, 
and if our sample was progressively contaminated by mis- 
classificd radio-loud quasars with redshift, it could account 
for the observed increase of Lg5 U(lm with redshift. Given the 
above considerations, this might happen because, for exam- 
ple, some of our high-redshift objects might be observed in 
a line of sight sufficiently close to the jet axis that the con- 
tribution from a more face-on quasar heated torus (which 
may be optically thick in the rest-frame mid-far infrared) 
becomes dominant at 850/im. 

One way to resolve this issue is to check directly whether 
unification works at submillimetre wavelengths, and we are 
currently undertaking a programme of SCUBA observations 
of radio-loud quasars to test this hypothesis. However, even 
if radio-loud quasars and radio galaxies are shown to differ 
at submillimetre wavelengths this need not invalidate the 
conclusions of this study, provided we can demonstrate that 
progressive quasar contamination with increasing redshift is 
not a problem. 

We can in fact test the likelihood of this selection ef- 
fect directly using the basic radio properties of the objects 
we have observed with SCUBA. In Figure ^, the rest-frame 
15 GHz core:lobe ratio is plotted against redshift and Lgso Mm 
for our sample. If the high-redshift objects in the sample are 
in fact buried radio-loud quasars masquerading as galaxies, 
they should be expected to display systematically higher ra- 
dio core:lobe ratios, owing to Doppler boosting of the radio 
core along the line of sight. The typical core:lobe ratio does 
not appear to change significantly with redshift, and so this 
analysis does not support the idea of a systematic bias to- 
wards buried quasars at high redshift within the SCUBA 
radio galaxy sample studied here. Moreover, there is no ev- 
idence that the more core dominated objects are more de- 
tectable with SCUBA. 

These mis-classified buried quasars may exist. Two of 
the objects in our sample, 6C0902+34 and 4C23.56, appear 
significantly more core-dominated than the rest, and are 
likely candidates since their core:lobe radio is comparable 
( ~ 10 per cent) to t hat typically displayed by radio quasars 
( |3pencer et al. 1991] ) . Furthermore, 6C0902+34 has a single- 
sided radio jet characteristic of a radio-loud quasar ( |Carilli| 



1995a ). In addition to the evidence at radio wavelengths, 



4C23.56 has a polarized u ltraviolet continuum consistent 
with scattered quasar light ( |Cimatti et al. 1998a ) . It is worth 
noting that the core:lobe ratio in these two sources is an or- 
der of magnitude higher than the typical value throughout 
the rest of the radio galaxy sample. Moreover, in any case, 
SCUBA failed to detect dust emission from either one. 
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Figure 8. Investigation of the rest-frame 15 GHz coredobe ratio 
for the objects in our sample. The top plot shows the core:lobe 
ratio against redshift. Solid circles represent the sources detected 
with SCUBA, and open circles represent the non-dctcctions. 
Down arrows indicate sources for which no radio core has been 
detected. In order to estimate the core flux at 15 GHz rest-frame, 
the core spectral index is required. For sources with more than 
one core detection, we have measured this spectral index. For 
sources with only one core detection, we have assumed the core 
to be flat. The bottom plot compares the coredobe ratio with the 
850-fim luminosity. An n o =1.0 cosmology is assumed; there is no 
significant difference between this plot and one for which f! o =0.1. 



7.3.2 Preferential selection of young radio sources at high 
redshift 

All realistic models of the time evolution of classical double 
radio sources (e.g. Kaiser et al. 1997; Blundell et al. 1999; 
Blundell & Rawlings 1999) predict that their lobe luminosi- 
ties decline with age once they are large enough to avoid syn- 
chrotron self-absorption. Older sources are thus (obviously) 
more likely to fall below the flux limit of a given survey if 
they reside at high redshift. A single flux-limited survey may 
thus be biased towards young radio sources at high redshift, 
whereas at lower redshifts a large range of radio source ages 
will be included in the same survey. 

The importance of this selection effect to our SCUBA 
survey is probably not a major concern for two reasons. 
First, the target sources are not selected from a single flux- 
limited sample, but rather are selected rather close to the 
flux-density limit of progressively deeper surveys - in other 
words at z ~ 4 the radio sources must be preferentially 
young, but the sources we have studied at, say, z — 1.5 
must also be fairly young to enter the brighter 3CRR sample 



Figure 9. The cumulative redshift distribution of the radio 
galaxy sample. The dotted line is for the entire sample, the solid 
line only considers the galaxies detected at 850 fim. For the de- 
tections, the median redshift of the distribution is < z >= 3.1. 




- CY+B99 minimum-z 
CY only minimum— z 
CY+B99 flat P(z) 
Borger et al. (1999) 



Figure 10. Figure 2 of Smail et al. 2000b. The full caption, 
as it appears in the original paper, is as follows: "Cumulative 
redshift distribution for the full submm sample. We have used 
the spectroscopic redshifts of those sources thought to be reliable 
(Table 1) and combined these with the probable redshift ranges of 
the remaining sources derived from their of 5 4 indices or limits. 
The solid line shows the cumulative distribution if we assume 
the minimum redshift distribution that is obtained if all sources 
lie at their lower z a limit given in Table 1 (dashed line is the 
equivalent analysis but restricted to the C Y models) . The effect of 
non-thermal radio emission, which drives down the Qj 5 ^ indices, 
means that this is a very conservative assumption if some fraction 
of the population harbor radio-loud AGNs. The dash-dotted line 
assumes a flat probability distribution for the sources within their 
z a ranges and a maximum redshift of z = 6 for those sources 
where we only have a lower limit on af 5 ® . Finally, the dotted line 
is the cumulative redshift distribution from Barger et al. (1999a) 
with two of the source identifications corrected as in Smail et al. 
(1999) and all blank- field/extremely red object candidates placed 
at z = 4." 



from which our lower-redshift targets are selected. This fact 
will undoubtedly weaken any correlation of source age with 
redshift (certainly at z > 1), compared with that expected 
from a single flux-limited survey as described by Blundell & 
Rawlings (1999). 

Second, even if some residual correlation between radio- 
source youthfulness and redshift remains in our sample, it 
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is not clear that this has any implications whatsoever for 
this submillimetre study. For it to be responsible for the 
trend of increasing Lgso M m with z, radio-source age must be 
intimately linked to starburst luminosity. Our data do not 
provide any evidence in support of this. For example, one of 
the largest radio sources, and hence possibly one of the old- 
est, 6C0905+39, was detected by SCUBA. To examine this 
more thoroughly, survival analysis was applied to test for 
correlations between T>n near (taken as indicative of the age 
of the radio source), z, and Lg5o Mm . As found for complete 
samples by Blundell et al. (1999), we find a negative cor- 
relation between D( inear and z. However, no evidence for a 
correlation between D;i near and Lgso^m was found (correla- 
tion tests yield a significance of ~ 40 per cent). We therefore 
find no internal evidence in support of the hypothesis that 
L850(im is closely linked to radio-source age, although the 
mildness of the Dn near '-z correlation means that caution is 
advised in attempting to in terpret Dunear as a reliable esti- 



mator of radio-source age dBlundcll et al. 199§ [Blundell &] 
RawlinlsTg^. 



8 CONCLUDING REMARKS 

In summary, our attempts to investigate and quantify the 
various possible biases which might conceivably afflict this 
study have simply served to reaffirm and strengthen our 
basic result, namely that the submillimetre luminosity of 
radio galaxies is primarily a function of redshift as illustrated 
in Figure (7J 

It therefore seems hard to avoid the straightforward 
conclusion that the observed increase in submillimetre de- 
tection rate and characteristic luminosity with redshift is 
due to the increasing youthfulness of the stellar populations 
of the radio galaxies in our sample. 

In a separate paper we will explore how the inferred evo- 
lution of gas mass and star-formation rate in these galaxies 
compares with the predictions of models of elliptical galaxy 
formation and evolution. However, it is interesting to briefly 
consider whether the apparently rather extreme evolution is 
peculiar to radio galaxies, or may in fact be typical of the 
cosmological evolution of dust and gas in massive ellipticals 
in general. 

In Figure ^| we show the cumulative redshift distribu- 
tion of our radio galaxy sample, along with the cumulative 
redshift distribution of the subset detected at submillime- 
tre wavelengths. This figure serves to re-emphasize that the 
high median redshift of our detected galaxies (z — 3) does 
not simply reflect the median redshift of the sample selected 
for observation (z — 2). However, what is particularly in- 
teresting is that the redshift distribution of our submillime- 
tre detected radio galaxies is statistically indistinguishable 
from current best estimates of the redshift distribution of 
sources detected in SCUBA surveys, as illustrated in Fig- 
ure [l[] ( |Smail et al. 2000b[ ). This comparison provides at 
least circumstantial evidence that the submillimetre evolu- 
tion of radio galaxies found here may indeed be symptomatic 
of the evolution of massive elliptical galaxies in general. 

At first sight, it may appear contradictory that the av- 
erage submillimetre luminosity of the radio galaxy sample 
continues to rise beyond z ~ 3 while the median redshift of 
the sample is z ~ 3. This results from the median redshift 



of our most luminous detections being higher at z ~ 3.6 for 
sources brighter than 5 mjy. 

This raises the interesting possibility that the most 
massive dust-enshrouded starbursts are confined to z > 3. 
Therefore the median redshift of submillimetre sources de- 
tected in blank-field surveys may prove to be a function of 
flux density. 

It will be some time before the redshift information 
in bright submillimetre surveys approaches that currently 
available for radio-selected samples. However, it will un- 
doubtedly be very interesting to see how this comparison 
evolves as submillimetre-selected samples are studied and 
refined in the years to come. 
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APPENDIX A: RADIO-SUBMILLIMETRE 
SPECTRAL ENERGY DISTRIBUTIONS 

We applied two fits to the radio spectra of our sample: 

1) Linear Fit - In log-space, a straight line is fitted to 
radio data with v > 1 GHz. For all galaxies in the sample, at 
least two data points satisfy this condition. Except for the 
noticeably curved spectra of 3C324, MG1744+18, 4C23.56, 
and MG1016+058, the straight-line fit seems to mimic the 
existing high-frequency data very well. 

2) Parabolic Fit - A physical model of synchrotron 
ageing could be applied to the spectra in order to determine 
how they might steepen. However, the fit is only required 
to accurately predict the shape of the expected steepening; 
the physics of the ageing synchrotron population do es no t 
need to be extracted from the fit. Blundell et al. ( 199£ ) 
have studied the properties and evolution of over 300 FRII 
radio sources selected from the 3C, 6C, and 7C surveys. 
They applied a Bayesian polynomial regression analysis to 
the fitting the flux densities. The process involved fitting 
polynomials of different degrees and associating a likelihood 
or probability to each. In some cases a first-order polyno- 
mial, or straight-line was the preferred choice. However, in 
most cases a second-order polynomial, or parabola, accu- 
rately modelled the curvature of the spectrum. A parabolic 
fit seems to be the best way of parameterising the curva- 
ture, as higher-order polynomials were never chosen by the 
regression analysis. Thus a polynomial of degree m — 2 was 
fitted, in log space, to the entire radio spectrum. 

T he b est-fit parameters of the two models are shown in 
Table 



Al 



For each galaxy in the sample, the radio-submillimetre 
spectral energy distribution (SED) is presented in Figure |XJ 
The SEDs are arranged in redshift order. The solid circles 
denote integrated continuum flux densities. The open cir- 
cles indicate radio core flux densities. For 3C356, two ra- 
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Table Al. Best-fit parameters (including the reduced \ 2 of the fit) of both straight-line and parabolic fits to the radio spectra (the fits 
were made in log-space). Note, the linear fits were applied to v > 1 GHz observations only. In addition, if only two data points were 
available (i.e. zero degrees of freedom), the slope and y-intercept of the line joining the two data points were calculated instead of using 
a fitting algorithm. No value of X 2 red * s gi ven m this case. Finally, all best-fit parameters are quoted to two decimal places; therefore 
errors less than 0.005 are quoted as being 0.00. 



Linear Fit: y = mx + k Parabolic Fit: y = ax 2 + bx + c 
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dio core identifications exist. The fainter core is denoted 
by open stars. If errors were unavailable, the error was 
taken to be 10 per cent of the flux density. Upper lim- 
its have been taken at the 3-a level if S/N< 3. The two 
models fit to the data are also displayed in Figure [X| 
The solid line is the parabolic fit to the spectrum, and 
the dash-dot line is the linear fit. The references for the 
data are: 3C277.2 - BRR00, LP; 3C340 - BRR00, LP; 
3C265 - BRR00, LP; 3C217 - BRR00, LP; 3C356 - 



BLR97, F93, LP; 3C368 - BLR97, BRR00, LP; 3C267 - 
BLR97, BRR00, LP; 3C324 - BCG98, BRB98, LP; 3C266 
- BLR97, BRR00, LP, LPR92; 4C13.66 - BWE91, GC91, 
BRR00, LML81, LRL, P90, WB92; 3C437 - BRR00, LP; 
3C241 - ASZ91, LP, vBF92; 6C0919+38 - 6C-II, CC98, 
FGT85, GC91, GSDC96; 3C470 - BLR97, LP; 3C322 - 
BRR00, LP; 6C1204+37 - 6C-II, CC98, DBB96, GC91, 
GSDC96, LLA95; 3C239 - BRR00, LP; 3C294 - BRR00, 
LP; 6C0820+36 - 6C-VI, CC98, DBB96, GC91, GSDC96, 
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LLA95; 6C0905+39 - 6C-II, CC98, DBB96, FGT85, 
GC91, GSDC96, LGE95; 6C0 901+ 35 - 6C-II, CC98, 
DBB96, GC91, GSDC96, NAR92(|l99§); 5C7.269 - BRROO, 
CC98, DBB96, GSDC96; 4C40.36 - C&C69, CMvB96, 
CRv097, FGT85, TMW79, WENSS; MG1744+18 - 
BWE91, CRv097, GC91, GSW67, KB17, LML81, P90, 
WB92; MG1248+11 - DBB96, LML81, LBH86, WB92; 
4C48.48 - 6C-V, CMvB96, CRv097, GSW67, TMW79, 
WENSS; 53W002 - 6C-II, HDR97, WBM91, WvHK84, 
OKS87, OvL87; 6C0930+38 - 6C-II, CC98, DBB96, 
FGT85, GC91, GSDC96, NAR92( |l992| ); 6C1113+34 - 6C- 
II, CC98, DBB96, GC91, GSDC96, LLA95A; MG2305+03 

- BWE91, GC91, GWBE95, LML81, P90, WB92; 3C257 - 
BWE91, CC98, DBB96, GC91, GSW67, GWBE95, HDR97, 
LML81, P90, vBS98, WB92; 4C23.56 - CMvB96, CRv097, 
P90, PS65, TMW79; 8C1039+68 - HWR95, Lthesis; 
MG1016+058 - BWE91, CC98, CFR98, DBB96, DSD95, 
GC91, GSW67, GWBE95, LML81, P90, WB92; 4C24.28 - 
CMvB96, CRv097, P90, PS65, TMW79, WYR96; 4C28.58 

- CC98, CMvB96, DBB96, PS65, TMW79; 6C1232+39 

- 6C-II, CC98, C&. K94 CRv097, DBB96, FGT85, GC91, 
GSDC96, NAR92( [l99^ ); 6C1159+36 - 6C-II, CC98, 
DBB96, GC91, GSDC96, LLA95; 6C0902+34 - 6C-II, C95, 
CC98, C&K94, COH94, DBB96, DSS96, GC91, GSDC96, 
HDR97, LLA95; TX1243+036 - CC98, CFR98, DBB96, 
GSW67, LML81, P90, vOR96, WB92; MG2141+192 - 
BWE91, CC98, CRv097, DBB96, GC91, GSW67, HDR97, 
WB92; 6C0032+412 - 6C-VI, BRE98, DBB96, FGT85, 
GC91, GSDC96, HDR97; 4C60.07 - 6C-V, CC98, CMvB96, 
CRv097, DBB96, TMW79, WENSS; 4C41.17 - 6C-VI, 
BC099, CC98, C&K94, CMvB90, COH94, DBB96, DHR94, 
FGT85, GSW67, TMW79, WENSS; 8C1435+635 - 6C-III, 
CC98, CRv097, DBB96, HDR97, HWR95, 195, IDHA98, 
LMR94; 6C0140+326 - 6C-VI, BRE98, CC98, DBB96, 
FGT85, RLB96. 
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Figure Al. Spectral Energy Distributions 
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Figure Al - continued 
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Figure Al - continued 
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